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Abstract 


The  eddy-correlation  measurement  technique  was  used  to  obtain 
micrometeorological  field  measurements  of  the  fluxes  of  momentum, 
sensible  heat,  and  ozone  at  3  and  8  meters  at  the  Pawnee  site.  A 
method  to  accurately  measure  mass  (ozone)  gradients  from  surface- 
layer-based  meteorological  towers  was  developed  and  used.  Both 
flux  and  gradient  measurements  were  used  for  the  determination  of 
eddy  diffusivities.  Ozone  was  used  as  the  mass  tracer  to  explore 
similar  relationships  between  the  eddy  diffusivities  of  momentum 
K^,  sensible  heat  Kh,  and  mass  Kc. 

Eddy-diffusivity  ratios  were  computed  using  dimensionless- 
gradient  ratios  classified  from  the  data  and  regression  models.  These 
ratios  were  classified  by  atmospheric  stability  determined  at  the 
geometric  mean  of  the  measurement  heights. 

The  assumption  of  similarity  between  the  eddy  diffusivities  of 
ozone  and  sensible  heat,  Kc  =  Kh,  based  on  scalar  turbulent  transfer 
theory,  was  verified  for  unstable  atmospheric  conditions.  Eddy 
diffusivities  of  sensible  heat  are,  however,  50%  greater  than  diffu- 
sivities of  ozone  for  stable  atmospheric  conditions.  Chemical  reac- 
tion of  ozone,  and/or  the  need  for  flux-measurement  corrections, 
decrease  the  resulting  values  for  ozone  diffusivities  during  stable 
periods.  Established  eddy-diffusivity  ratios  for  water  vapor  and 
momentum  are  valid  for  ozone  and  momentum  under  stable  atmo- 
spheric conditions  over  smooth  terrain,  but  not  under  unstable 
conditions  for  flow  disturbed  by  irregular  terrain. 

The  relationships  between  the  eddy  diffusivities  of  momentum 
and  ozone,  as  well  as  those  between  momentum  and  sensible  heat, 
are  controlled  by  free-convection  conditions,  <  Kc  and  <  K^; 
these  results  are  inconclusive  for  unstable-atmospheric,  smooth- 
terrain  conditions  for  the  Pawnee  site  during  easterly  winds.  For 
disturbed-flow  conditions  at  the  Pawnee  site  during  westerly  winds, 
the  eddy  diffusivities  for  momentum  were  enhanced  and  both  K^/K^ 
and  Kh/Km  become  approximately  constant  (0.35  and  0.5,  respec- 
tively) for  both  stable  and  unstable  conditions. 

Theoretical  flux  measurement  corrections  were  applied  to  a 
subset  of  the  data  in  an  attempt  to  investigate  the  concept  of  a 
constant  flux  layer  and  to  validate  flux  correction  theory.  The  results 
indicate  that  flux  corrections  are  warranted.  However,  improve- 
ments to  flux  correction  theory  are  needed,  especially  for  stable 
atmospheric  conditions,  as  demonstrated  by  the  significant  increase 
in  data  variability  as  a  result  of  these  corrections.  The  "constant  flux- 
layer"  concept  holds  for  momentum  flux  and  sensible  heat  flux,  but 
is  questionable  for  ozone  flux. 
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Eddy  Diffusivities  for  Sensible  Heat,  Ozone,  and  Momentum 
From  Eddy  Correlation  and  Gradient  Measurements 

K.  F.  Zeller 


MANAGEMENT  IMPLICATIONS 

Fluxes  of  trace  gases  ("greenhouse  gases")  in  the 
atmospheric  boundary  layer  have  attracted  the  con- 
cerned interest  of  the  scientific  community  because  of 
the  major  climate  changes  projected  for  the  present  and 
the  near  future.  Although  concentrations  of  these  trace 
gases  are  known  and  have  been  monitored  for  years,  it 
is  information  on  sources,  sinks,  and  transport  of  these 
gases  that  is  necessary  for  predicting  future  concentra- 
tions and  the  consequence  of  climate  change.  The  mea- 
surement and  analysis  of  fluxes  and  vertical  gradients  of 
ozone,  as  well  as  traditional  micrometeorological  fluxes 
of  heat  and  momentum  over  long  time  periods,  offer 
comparisons  and  better  understanding  of  measured 
boundary  layer  eddy  diffusivities  of  momentum,  heat, 
and  ozone.  These  measurements  over  dry  grasslands 
provide  region-specific  fluxes,  information  on  the  trans- 
port behavior  of  ozone,  and  the  documentation  of  how 
to  make  such  flux  measurements. 

INTRODUCTION 

Since  the  mid  1960's  with  the  advent  of  modern- 
environmental  legislation,  the  horizontal  transport  and 
dispersion  characteristics  of  the  wind  in  the  atmo- 
spheric surface  layer  have  been  studied  and  modeled 
with  the  goals  not  just  of  basic  research,  but  of  environ- 
mental management  and  protection.  Modern  environ- 
mental decisions  have  been  based  on  modeled  and 
measured  pollutant  concentrations  within  the  ambient 
air  (Cermak  1975;  Kao  1984;  Randerson  1984).  Now,  as 
a  result  of  the  "acid  rain"  issue,  emphasis  is  being 
shifted  to  the  amount  of  a  pollutant  delivered  to  and 
"taken  up"  by  various  components  of  the  earth's  sur- 
face: bare  soil,  water  surfaces,  biological  canopy,  or 
man-made  or  natural  material  surfaces  (NRC  1983).  The 
amount  of  a  pollutant  removed  to  a  surface — mass  per 
unit  area  per  unit  time — is  the  pollutant  deposition  flux. 
Knowledge  of  and  predictions  of  pollutant  deposition 
are  required  for  acid  rain  related  and  other  air  pollutant 
related  impact  assessments.  Deposition  has  consequently 
become  an  important  consideration  in  recent  modeling 
developments  (Hicks  1984;  NCAR  1985;  Venkatram  et 
al.  1988;  Walcek  et  al.  1986;  Walcek  and  Chang  1987). 
Deposition  can  be  either  "wet"  or  "dry,"  depending 
upon  the  delivery  process.  An  air  pollutant  making  its 
way  to  a  surface  without  the  aid  of  precipitation  would 
be  dry-deposited  upon  leaving  the  atmosphere. 


The  physical  origin  of  the  deposition  of  a  pollutant, 
c,  can  be  understood  by  visualizing  its  mass  balance 
within  any  airshed  as  a  simple  control  volume  (McRea 
and  Russell  1984): 

—  I  cdv  +    f  Fr  .  nds  -    |  Rdv  =  0  [1] 

*  Jv         Js  Jv 

where  the  first  term  accounts  for  emissions  of  c  [gm  m"3 
or  parts  per  billion  (ppb)  by  volume]  integrated  within 
the  control  volume,  V;  the  second  term  accounts  for  the 
flux  Fc  [gm  m"2  s"1] ,  of  c  integrated  across  the  boundaries 

S  (  n  is  the  unit  normal)  of  the  surface;  and  R  is  the  rate 
of  chemical  production  or  reduction  integrated  within 
the  control  volume. 

Assuming  that  the  airborne  pollutant  in  question  is 
nonreactive  and  that  it  eventually  exits  the  atmosphere 
at  the  earth's  surface,  the  correct  evaluation  of  the 
surface  integral  (second  term)  in  [ll  over  the  lower 
control  volume  boundary  determines  the  ultimate  cor- 
rectness of  the  modeled  scenario.  The  flux  of  material 
for  a  specified  pollutant  species  c  passing  through  the 
lower  surface  is  its  deposition  flux  Fc. 

Governing  Equations 

The  expression  for  the  concentration  of  a  pollutant 
(scalar  quantity)  transported  through  a  boundary  layer 
is  (Businger  1986;  Rohsenow  and  Choi  1961): 

9t      axj  5x2 

1 

where 

c  =    c(t,  Xj),  the  instantaneous  concentration  (in 
gm  m"3  or  ppb)  (Einstein's  summation  con- 
vention applies  over  repeated  indices); 
Uj  =   u(t,  Xj),  the  instantaneous  orthogonal  compo- 
nents of  wind  velocity; 

Dc  =  the  molecular  diffusivity  of  c;  and 

S  =   a  source  or  sink  of  c. 

When  the  instantaneous  concentration  and  instan- 
taneous velocities  are  decomposed  into  their  average, 
c ,  and  fluctuating  or  turbulent  parts,  c',  (i.e.,  c  =  c  +  c' 
and  Uj  =  ui  +  u\),  and  then  combined  and  averaged,  [2] 
becomes: 
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[3] 


Assuming  steady-state,  horizontal  homogeneity,  and  no 
local  sources  or  sinks  allows  the  first,  second,  and  fifth 
terms  of  [3]  to  be  canceled  leaving: 


3w'c' 
3z 


=  D 


3^c 
C3z2 


[4] 


where  w'  =  u'3  and  z  =  x3,  the  vertical  wind  velocity 
fluctuation  and  vertical  distance  coordinate,  respec- 
tively. Integration  of  [4]  between  the  earth's  surface  and 
height  z  yields  a  description  of  the  vertical  flux  of  c 
between  those  levels: 


w'c 


w  c 


=  D, 


'acp 

3z 


D, 


'acp 

3z 


[5] 


'o 


Since  w'  approaches  zero  at  the  surface,  (w'c')o  is 
minimal  and  can  be  disregarded.  Also  (Wc')z,  the  verti- 
cal turbulent  transport  at  height  z  is  many  orders  of 
magnitude  greater  than  Dc(3c  /  3z)z.the  vertical  mo- 
lecular diffusion  transport  at  height  z;  therefore,  it  can 
also  be  disregarded.  This  leaves  the  expression: 

=  -bJ— )  [6] 
A) 


3z 


to  describe  the  vertical  transport,  flux,  or  dry  deposition 
of  pollutant  c  (Businger  1986).  Equation  [6]  is  valid  for 
any  scalar  (e.g.,  temperature,  water  vapor)  given  the 
same  assumptions  in  the  derivation.  The  surface  mo- 
lecular transport,  second  term  in  [6],  cannot  be  mea- 
sured directly  (Businger  1986);  therefore,  the  first  term, 
wV,  measured  at  height  z,  is  the  only  opportunity  to 
directly  measure  Fc  ,  the  vertical  flux  of  c. 


Fc  =  w'c' 


[7] 


The  eddy  correlation  technique  is  used  to  measure  the 
right-hand  side  of  [7]  directly,  although  indirect  mea- 
surement techniques  exist. 


Eddy  Diffusivities 

Specification  of  flux  for  momentum  t,  mass  Fc, 
sensible  heat  H,  and  water  vapor  XE  have  historically 
been  based  on  the  gradient  (or  profile)  approach  (Bache 
1986;  Munn  1966;  Rosenberg  1974;  Sellers  1965).  The 
assumption  is  that  the  same  type  of  equation  for  molecu- 
lar diffusion  transfer  applies  to  turbulent  transfer: 


t=  PKm(3u73z) 
Fc=  -Kc  (3c/3z_) 
H=  -pCpKhO0/3z) 
XE  =  -pLKv  (3v73z) 


[8a] 
[8b] 
[8c] 
[8d] 


where  Km  is  the  eddy  (or  turbulent)  diffusion  coeffi- 
cient for  momentum,  Kc  for  ozone  (mass),  Kh  for  sen- 
sible heat,  and  for  water  vapor;  p  is  the  density  of  air; 
v  is  the  average  water  vapor;  and  6  is  the  average 
potential  temperature.  The  concept  of  an  eddy  diffusion 
coefficient  for  momentum  is  rooted  in  Prandtl's  "mixing 
length  theory"  (Haltner  and  Martin  1957). 

Two  recent  models  have  used  the  K  approach.  The 
Regional  Acid  Deposition  Model,  RADM  (Chang  et  al. 
1986;  NCAR  1985),  developed  for  EPA  at  NCAR,  for 
instance,  integrates  [8b]  as  a  starting  point  in  the  devel- 
opment of  a  method  to  calculate  pollutant  fluxes.  In 
RADM,  Kc  is  provided  for  through  a  similarity  relation- 
ship with  Km.  The  Acid  Deposition  and  Oxidant  Model, 
ADOM  (Venkatram  et  al.  1988),  an  internationally  spon- 
sored model  similar  to  RADM,  uses  the  same  approach. 

The  recommended  engineering  approach  to  speci- 
fying Kc  for  ozone  or  for  any  other  airborne  mass  is  to 
assume  the  principle  of  similarity  (Galbally  1971;  Munn 
1966;  Oke  1978;  Rider  and  Robinson  1951;  Sellers  1965): 


Kc  =  Kh  =  Kv  =  aKm 


[9] 


where  a  is  a  constant  (~  1.35)  for  stable  atmospheric 
conditions  and  a  function  of  stability  for  unstable  atmo- 
spheric conditions  (Businger  1986). 


Goals  and  Hypotheses 

This  paper  investigates  exchange  processes  between 
the  atmosphere's  surface  layer  and  the  earth's  surface. 
Pollutants  are  transported  from  the  earth's  air-surface 
layer  to  an  underlying  surface.  Measurements  of  sur- 
face-layer fluxes  and  vertical  gradients  of  sensible  heat, 
ozone,  and  momentum  were  acquired  over  a  relatively 
flat-terrain,  shortgrass  steppe  prairie  (Zeller  et  al.  1989). 
Diffusivities  obtained  directly  from  these  measurements 
are  compared  to  established  empirical  relationships, 
and  the  similarity  assumptions,  [9],  are  evaluated. 

Mass  diffusivities,  such  as  the  diffusivity  for  ozone, 
are  used  by  the  scientific  and  engineering  communities 
to  describe  the  diffusion  of  mass  through  the  atmo- 
sphere. The  results  presented  in  this  paper  can  be  used 
to  improve  pollution  impact  assessments  and  to  help 
establish  relationships  between  pollutants  and  the  dam- 
age they  cause.  More  clearly  understood  and  specifi- 
cally defined  mass  diffusivities,  for  example,  will  more 
accurately  account  for  actual  mass  deposition  of  pollut- 
ants. 

This  study  had  two  main  objectives: 

1.  To  investigate  the  validity  of  using  eddy  diffu- 
sivities of  sensible  heat  or  momentum  to  specify 
the  eddy  diffusivity  for  ozone  under  various 
atmospheric  stabilities. 


2 


2.  To  investigate  the  assumption  of  a  constant 
vertical  flux  layer  for  heat,  momentum,  and 
ozone  within  the  lower  atmospheric  surface 
layer. 

My  aim  was  to  investigate  time-honored  assump- 
tions that  are  commonly  used,  but  as  yet  have  not  been 
completely  verified  by  experimental  methods.  In  fact, 
attempts  to  prove  or  disprove  these  assumptions  have 
met  with  mixed  results  (Businger  1986).  The  data  col- 
lected at  the  Pawnee  site  were  culled  and  edited  to 
provide  representative  data  records,  which  were  then 
used  to  calculate  and  analyze  dimensionless  gradient, 
eddy  diffusivity,  and  atmospheric  stability  quantities. 
Objective  1  was  evaluated  by  comparing  measured  re- 
sults and  regression  models  of  the  measured  results 
with  accepted  empirical  theories  for  smooth  terrain  at 
various  atmospheric  stabilities.  Objective  2  was  evalu- 
ated by  comparing  pairs  of  observed  corrected  and 
uncorrected  fluxes  at  two  levels  with  the  expected  ratio 
of  1:1. 


SUMMARY  OF  PREVIOUS  INVESTIGATIONS 
Eddy  Diffusivities 

Flux  equations  [8a-d]  are  the  turbulent  analogies  to 
the  basic  physical  law  relating  molecular  viscosity  to 
molecular  diffusion.  The  eddy  diffusivity  (also  called 
eddy  viscosity,  eddy  exchange  coefficient,  and  austausch 
coefficient)  Km  is  the  turbulent  analogy  to  Newton's  law 
for  molecular  viscosity.  The  eddy  diffusivities  Kc,  Kh, 
and  are  the  eddy  diffusivities  analogous  to  Fourier's 
and  Fick's  laws  for  the  molecular  diffusion  coefficients 
for  mass,  heat,  and  water  vapor.  The  units  of  these 
coefficients  are  the  same,  cm2s_1;  however,  the  eddy 
diffusivities  are  several  orders  of  magnitude  greater 
than  their  molecular  counterparts.  Eddy  diffusivities 
are  based  on  assumptions  about  the  similarity  between 
molecular  and  turbulent  transfer,  not  on  sound  physical 
laws  (Arya  1988).  They  are  properties  of  the  flow,  not 
simply  fluid  properties.  Therefore,  variations  in  mea- 
sured Kq's  are  expected  as  the  atmospheric-surface- 
layer  characteristics  change. 

Based  on  scaling  arguments,  Lettau  (1951)  limits  the 
value  of  K  to  a  maximum  of  108  cm2s_1,  and  he  predicts 
that  Km  will  vary  from  ~1  to  -10  cm  s  between  the 
earth's  surface  and  1-km  height.  Jacobi  and  Andre 
(1963)  provide  estimated  height-dependent  values 
(cm2s"a)  based  on  atmospheric  stability: 


Height 


3  m 
8  m 


Inversion 
30 
80 


Stability 

Stable  Neutral 
250  5,500 
790  16,000 


Unstable 
55,000 
160,000 


Given  temperature  and  mass  are  both  scalars,  heat 
and  mass  transfer  by  the  turbulent  atmosphere  should 
be  similar,  hence  the  expected  ratio  Kh/Kc  is  1.  Empiri- 
cal data  for  various  flows  indicate  that  the  ratio  Kh/1^= 
a  should  be  in  the  range  0.5  and  1  (Bird  et  al.  1960); 
however,  based  on  atmospheric  experiments,  a  de- 
pends on  stability  and,  as  presented  above,  is  sometimes 
estimated  as  1.35  for  stable  conditions.  Businger  et  al. 
(1971)  provide  an  empirical  relationship  for  the  ratio  a 
as  a  function  of  atmospheric  stability: 

1.35(1  -9Q1/2    .      r    n  ,  v 
for  £<0  (unstable) 


a  = 


a  = 


(1-15QJ 
(l  +  4.7z) 


[10] 


for  z  >  0  (stable ) 


(0.74  +  4.7z) 
where 

£  =    z/L,  an  atmospheric  stability  parameter; 

z  =    height  above  ground;  and 

L  =   the  Monin-Obukhov  length,  [13]. 

Based  on  [10],  a  is  equal  to  1.35  only  when  £— 0, 
then  it  quickly  drops  to  -1  for  £>0.  An  earlier  empirical 
result  given  by  Swinbank  (1968),  a=2.7£025,  for  £<0 
better  describes  the  Pawnee  site  results. 

In  an  attempt  to  preview  the  behavior  of  Kc/Km, 
Y^/Km  is  used  here  to  represent  a  "mass"  diffusivity  as 
opposed  to  Kh/Km.  Rosenberg  (1974)  provides  an  em- 
pirical relationship  for  the  ratio  of  Kv/Km  as  a  function 
of  atmospheric  stability  based  on  experiments  by  Pruitt 
et  al.  (1971): 


[U] 


Ky/I^  =  1.13(l-60Ri)0  074  Ri<0.017  (unstable) 

fc^/^  =  1.13(l+95Ri)-011  Ri<-0.011  (stable) 
where 

Ri  =  the  gradient  Richardson  number, 
another  stability  parameter,  [14]. 

Note  that  an  asymptotic  overlap  exists  in  [11]  around 
Ri=0.  Panofsky  and  Dutton  (1984)  use  the  approxima- 
tions £=Ri  for  £<0  and  £=Ri/(l-5Ri)  for  £>0.  Equations 
[10]  and  [11],  the  ratio  Kh/Kv,  [10]  divided  by  [11],  and 
a=2.75£  are  plotted  in  figure  1  as  functions  of  £. 
Figure  1  shows  that  the  eddy-diffusivity  similarity  as- 
sumption, [9],  does  not  hold  exactly  based  on  data 
reported  in  the  literature,  as  the  empirically  determined 
ratio  K^/Ky,  is  not  precisely  1.  The  ratio  K^/K^  reaches 
a  constant  value  of  1.2  for  £>1  and  is  still  increasing  at 
the  value  1.8  for  £<-5.  The  ratio  K^/fC^  plotted  in  figure 
1 ,  different  from  unity,  is  a  good  preview  estimate  based 
on  reported  data  available  for  Kh/Kc  because  both  Kc 
and  Ky  are  mass  diffusivities  different  from  and  Kh, 
momentum  and  heat  diffusivities. 

Eddy-diffusivity  similarity,  [9],  is  used  throughout 
the  technical  literature  (e.g.,  Businger  1986).  A  variety 
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Figure  1 — Empirically  derived  eddy  diffusivities  vs.  stability 
parameter  (z/L). 

of  techniques  can  be  used  to  estimate  Kc.  For  instance, 
by  specifying  Kh  =  K^,  H  and  X.E  can  be  determined 
through  [8c]  and  [8d]  using  the  Bowen  ratio  method 
(Rosenberg  1974),  which  involves  measuring  the  sur- 
face energy-balance  components  using  well-proven, 
reliable  instruments.  Once  XE  or  H  are  known,  either 
[8c]  or  [8d]  with  either  calculated  Kh  or  calculated 
can  be  used  to  specify  Kc  by  [9]  (Sellers  1965).  In 
another  approach,  X.E  and  H  are  determined  by 
calculating  t  using  the  logarithmic  wind  profile  and 
assuming  the  similarity  relationship  with  aKm  and  the 
other  eddy  diffusivities,  [9]. 

The  mass  eddy  diffusivity  Kc  and  flux  Fc  can  also  be 
estimated  by  applying  a  modified  Bowen  ratio,  a  gradi- 
ent technique  that  involves  measuring  dc/dz,  dd/dz, 

Table  1 .  —  Ozone 


and  H,  then  using  [8b]  with  the  similarity  assumption. 
This  approach  has  been  used  by  many  (Droppo  1985; 
Woodruff  1986);  however,  concentration  differences 
measured  over  short  vertical  distances  are  difficult  to 
resolve  because  of  limited  instrument  resolution.  To 
circumvent  this  limitation,  one  instrument  is  typically 
used  with  the  sampling-intake  height  switched  periodi- 
cally between  the  chosen  8z  levels  over  the  data  averag- 
ing period  (Businger  1986).  Woodruff  (1986)  quantified 
errors  up  to  32%  using  this  approach  to  the  gradient 
technique. 

Table  1,  a  review  of  experimental  results  from  cur- 
rent literature,  shows  a  paucity  of  dc/dz  data  are  avail- 
able from  which  to  estimate  Kc.  This  paper  includes  the 
development  of  a  better  technique  to  measure  dc/dz 
specifically  for  ozone  that  will  be  applicable  to  other 
trace  gases. 

All  of  the  methods  that  employ  [8b]  for  determining 
Kc  and  Fc  use  the  similarity  assumption  [9].  All  of  the 
acid  deposition  models  developed  to  date  use  the  simi- 
larity assumption,  although  proof  of  this  assumption 
through  experiment  has  not  been  thoroughly  explored. 


Eddy  Correlation  Flux  Measurements 


The  direct  measurement  of  the  flux  w'c',  [7],  is 
called  the  eddy  correlation  technique.  It  has  long  been 
considered  the  most  reliable  (Arya  1988)  and  most 
fundamental  approach  to  flux  measurements  as  mass, 
heat,  and  momentum  are  primarily  transported  by  tur- 
flux  experiments. 


Reference 

Method* 

Sample 
height 
location 

Instruments 

Sample 
freq. 
(Hz) 

EC: 
(ppb  m  s1) 

Measured  range 

GR:  Vd: 
(ppb  rrr1)       (cm  s"1) 

Delany  et  al. 
(1986) 

EC 

4  &  8  m 

Chemiluminescence, 
sonic  anemometer 

10 

0  to  -0.7 

Droppo  (1985) 

EC 
GR 

9.1  m 
0.75  &  5.5  m 

Chemiluminescence, 
uvw  Gill  anemometer 

? 

0  to  -0.55 

-0.55  to  3.38 

Lenschow  (1982) 

EC/Vd 

Aircraft 

Chemiluminescence 

? 

-0.05;  1.0 

Meyers 

&  Yuan  (1987) 

EC/Vd 

? 

Summary  report 

? 

0  to  -0.04 
(daily  average) 

Wesely  et  al. 
(1982) 

EC/Vd 

5.2  m 

Chemiluminescence, 
assume:  uvw  Gill 

20 

-0.05  to  -0.3 

Wesely  et  al. 
(1981) 

EC 

5  to  8  m 
(winter) 

Chemiluminescence, 
uvw  Gill 

? 

-0.01  to  -0.13 

Wesely  et  al. 
(1978) 

EC 

4  &  5  m 

Chemiluminescence, 
uvw  Gill 

? 

-0.2  to  -0.55 

Zeller  et  al. 
(1989) 

EC 

6  m 
uvw  Gill 

Chemiluminescence, 

14 

0  to  -0.5 

0  to  -0.7 

Zeller  &  Hazlett 
(1989) 

EC 

6  m  Chemiluminescence, 
sonic  anemometer 

14 

0  to  -0.3 

0  to  -0.35 

GR:  gradient 

EC:    eddy  correlation  flux 


ECA/d:     eddy  correlation  with  only  deposition  velocity  reported  (Vd=Fc/  C ) 


not  specified  in  report. 


4 


bulent  motion.  The  eddy  correlation  method  also  pro- 
vides the  means  to  measure  sensible  heat  flux  H,  and 
momentum  flux  r,  directly: 

H  =  pCp  w7^7  [12a] 

t  =  p  u'w'  [12b] 

Requirements  for  valid  eddy  correlation  measure- 
ments are  stringent.  Until  the  recent  advent  of  fast- 
response  chemical  sensors  and  real-time  data  process- 
ing, eddy  correlation  was  a  preferred  method  in  theory 
only.  The  data  sampling  rate  has  to  be  fast  enough  to 
measure  most  of  the  eddies  contributing  to  the  flux. 
According  to  Kaimal  (1975)  and  his  co-workers  (Kaimal 
et  al.  1972)  a  sampling  frequency  of  about  10  Hz  or 
greater  is  needed  for  typical  surface-layer  measure- 
ments that  they  consider  valid  for  sampling  heights 
above  3  meters.  Below  3  meters  the  larger  transporting 
turbulent  eddies  are  restricted  by  the  earth's  surface, 
and  some  degradation  in  the  measured  flux  can  be 
expected  because  of  instrument  resolution  when  sonic 
anemometers  are  used  (Haugen  et  al.  1971;  Kaimal 
1969).  Although  in  theory  ensemble  averages  are  as- 
sumed when  developing  working  [3]-[7],  in  practice 
they  cannot  be  achieved  and  are  replaced  by  time 
averages  assuming  homogeneous  and  stationary  turbu- 
lence (Panofsky  and  Dutton  1984).  Based  on  field  expe- 
rience, the  averaging  period  for  eddy  correlation  data 
collection  should  be  between  0.25  and  1  hour;  short 
enough  to  maintain  stationarity,  but  long  enough  to 
capture  the  dominant  energy  contributing  eddies  that 
contribute  to  the  surface-layer  flux  (Brutsaert  1988). 
Sampling  details  are  examined  in  greater  detail  in  ap- 
pendix A. 

To  sample  at  10  Hz  or  faster,  fast-response  sensors 
with  short-time  responses  are  required  (Arya  1988). 
Restrictions  on  sensor  size,  placement,  and  sensing- 
element  geometry  also  exist  (Businger  1986).  In  cases 
where  sensors  are  not  exactly  matched,  are  not  collo- 
cated, do  not  sample  a  single  point,  or  cause  flow 
distortions,  flux  measurements  can  fall  short  of  true 
values.  In  such  cases,  theoretically  based  transfer  func- 
tions can  be  used  to  correct  flux  measurements  to  obtain 
the  correct  flux  (Kristensen  and  Fitzjarrald  1984;  Moore 
1986;  Zeller  et  al.  1989).  (These  transfer  coefficients  are 
detailed  in  appendix  A,  Zeller  et  al.  1989;  Zeller  1990). 
In  practice,  however,  correcting  flux  data  can  add  con- 
siderable scatter  to  the  results.  It  is  therefore  appropriate 
to  design  flux  measurements  in  such  a  way  as  to  elimi- 
nate these  potential  problems  as  much  as  possible 
(Wyngaard,  1989a). 

The  assumptions  that  lead  to  [4]  require  that  eddy 
correlation  measurements  be  made  above  flat,  uniform 
terrain.  In  practice,  flat,  expansive,  uniform  terrain  will 
provide  an  adequate  fetch  so  that  horizontal  homogene- 


ity can  be  assumed  within  a  reasonable  degree  of  error. 
Because  most  sites  are  not  perfectly  uniform  and  an 
ensemble  average  cannot  be  obtained,  the  vertical  com- 
ponent in  the  second  term  of  [3]  is  not  always  zero  (i.e., 
w  #0).  Recent  theoretical  developments  have  allowed 
for  some  relaxation  of  this  requirement  by  employing 
coordinate  rotation  forcing  w  =  0  (Wesely  1970),  thus 
adjusting  the  measured  eddy  correlation  to  the  value 
perpendicular  to  the  mean  calculated  wind-flow 
streamline,  (which  is  not  always  parallel  to  the  earth's 
surface)  for  the  sample  period.  This  approach  retains 
the  validity  of  [7]  and  has  been  found  reliable  for 
calculated  streamlines  within  ±5  degrees  of  the  surface 
plane  (McMillen  1986). 

Within  the  last  decade  several  researchers,  includ- 
ing this  writer,  have  attempted  to  measure  ozone  fluxes 
by  the  eddy  correlation  and  the  gradient  methods.  Table 
1  gives  a  representative  list  of  these,  including  the 
investigator,  method,  sampling  height,  instruments  used, 
sampling  frequency,  and  general  results  (order  of  mag- 
nitude) in  terms  of  either  Fc  for  flux,  dc/dz  for  mass 
gradient,  or  Vd  =  Fc/c  for  deposition  velocity.  Data 
reported  herein  are  within  the  same  orders  of  magnitude 
as  reported  in  previous  studies  (table  1). 

Surface-Layer  Considerations 

As  stated  previously,  eddy  diffusivities  describing 
turbulent  transfer  are  functions  of  the  fluid  flow.  Under- 
standing the  vertical  flux  or  dry  deposition  of  a  mass  to 
a  surface,  such  as  ozone  to  a  shortgrass  steppe  prairie 
canopy,  therefore,  requires  a  knowledge  of  the  flow 
properties  or,  in  this  case,  the  micrometeorological 
variables  just  above  the  surface  being  investigated. 
Monin-Obukhov  (M-O)  surface-layer  similarity  was  the 
micrometeorological  concept  used  to  provide  insight 
for  studying  the  behavior  of  surface-layer  turbulence 
(Wyngaard  1973).  The  basic  premise  for  mean  flow  and 
turbulent  characteristics  proposed  by  Monin  and 
Obukhov  in  1954  for  a  horizontally  homogeneous  sur- 
face is  that  they  depend  on  only  four  variables:  z,  the 
height  above  the  surface;  u*,  the  friction  velocity; 
H/pCp,  the  surface  kinematic  heat  flux;  and  g/T,  a 
buoyancy  variable.  Here  g  is  the  acceleration  due  to 

gravity  and  u*,  the  friction  velocity,  defined  as 
 1/2 

(-u'w')  •  As  these  four  M-O  variables  are  constituted 
by  only  three  fundamental  dimensions — length,  time, 
and  temperature — Buckingham's  PI  theorem  demands 
that  only  one  independent  dimensionless  parameter 
can  be  formed  by  them.  Obukhov  chose  the  buoyancy 
length  scale,  L,  in  1946  as  the  most  appropriate  (Arya 
1988): 

-u3T 

L  =  =  [13] 

kg  we' 
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where  w'0'  =  H/pCp.  Since  w'T'  was  measured  at  the 
Pawnee  site  w'0'  — 1.05  w'T'  for  1,641  meters  elevation 
is  used  to  calculate  L.  When  the  height  aboveground,  z, 
is  small  compared  to  L,  mechanical  turbulence  gener- 
ated by  wind  shear  prevails.  For  z  >  |L| ,  buoyancy 
effects  start  to  outweigh  mechanical  effects  (Clarke  et 
al.  1982).  Hence  z/L  =  £  is  an  important  buoyancy  or 
atmospheric  stability  parameter  and  is  traditionally 
used  in  most  modern  micrometeorological  studies.  As 
u*  and  H  have  been  difficult  to  measure  directly,  the 
gradient  Richardson  number  (Ri)  stability  parameter  is 
often  measured  and  used  to  indirectly  calculate  £  through 
the  empirical  relationships  given  in  Panofsky  and  Dutton 
(1984)  Chapter  II: 
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[14] 


The  M-0  similarity  theory  predicts  that  the  dimen- 
sionless  wind  shear,  the  dimensionless  sensible  heat 
gradient,  and  the  dimensionless  mass  gradient  are  all 
universal  functions  of  £  (Arya  1988;  Busch  1973): 


kz  3u 


=  U0 


kz  30  fn 
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[15a] 
[15b] 
[15c] 


Here  k  =  von  Karman  constant,  a  universal  constant 
empirically  derived  from  the  "law  of  the  wall"  when 
4>m(0)  =  1; 
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the  surface-layer 
temperature  scale; 

the  surface-layer  ozone 
concentration  scale; 

+  T     the  potential  temperature 
/    gradient;  and 

the  adiabatic  lapse  rate. 


The  empirical  forms  for  <J>m  and  <J>h,  determined 
from  many  field  studies,  are  the  subject  of  past  and 
ongoing  research  (Brutsaert  1988;  Fazu  and 
Schwerdtfeger  1989;  Hogstrom  1988;  Pruitt  et  al.  1971; 
Yaglom  1977): 


(1-b^)'174      for  £<0  (unstable) 


for  £>0  (stable) 


b3(l-b4n-1/2  for^<0 


[16a] 


[16b] 


b1+b2^ 


for  £>0 


The  relationships  established  from  the  1968  Kansas 
Experiment  are  generally  used  (Arya  1988;  Businger  et 
al.  1971)  and  will  be  used  here  to  compare  the  Pawnee 
site  results  with  existing  empirical  theory.  The  follow- 
ing values  were  established  for  a  flat  grassland  canopy 
and  represent  the  case  for  smooth  terrain:  b1  =  15; 
b2  =  4.7;  b3  =  0.74;  and  b4  =  9.  Typically,  <$>c  is  taken  to 
be  equal  to  <J>h  which  again  is  the  similarity  assumption. 

Velocity,  temperature,  and  mass  profiles  are  deter- 
mined by  integrating  [15a-c]  with  respect  to  height,  z 
(Arya  1988): 
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where  zQ,  the  surface  roughness,  is  a  constant  of  integra- 
tion evaluated  at  u  =  0.  It  is  typically  smaller  in  value 
compared  to  the  physical  height  of  the  surface  elements 
causing  it.  The  temperature  and  mass  concentration  at 
z0  are  0O  and  cQ.  The  functions  ^[Q  are  similarity 
functions  related  to  (J)^)  (Arya  1988): 
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The  above  relationships  for  u  and  i|jm  in  [17a]  and  [18b] 
are  used  to  establish  z0  for  the  Pawnee  study  site  dis- 
cussed in  this  paper. 


EXPERIMENTAL  DESIGN 
Field  Site 

The  data  analyzed  were  obtained  at  the  USDA 
Forest  Service  Pawnee  eddy-correlation  research  site. 
The  Pawnee  site  (fig.  2)  is  50  km  northeast  of  Fort 
Collins,  Colorado,  at  the  Central  Plains  Experimental 
Range  (CPER),  a  research  area  adjacent  to  and  west 
of  the  Pawnee  National  Grasslands.  The  CPER  is  a 
6,500-ha  area  of  shortgrass  prairie  that  is  administered 
by  the  USDA  Agricultural  Research  Service.  Also  lo- 
cated at  the  CPER  is  a  long-term  ecological  research 
(LTER)  station  managed  by  Colorado  State  University 
(CSU)  and  funded  by  the  National  Science  Foundation. 
The  exact  location  of  the  eddy-correlation  measurement 
site  is  within  a  6-ha  grazing  exclosure  in  Weld  County: 
section  27, 1,641  m  elevation,  latitude  40°28'23"  N,  and 
longitude  104°45'15"  W.  The  exclosure  was  established 
in  1969  as  an  appropriate  and  "representative"  site  for 
meteorological  measurements.  The  area  is  mostly  level, 
bordered  on  all  sides  by  gently  undulating  plains,  and  is 
not  near  major  ridges  or  bluffs  that  would  seriously 
influence  measurements.  There  is,  however,  a  drainage 
swale — Little  Owl  Creek — that  runs  north-south  ca.  400 
m  west  of  the  meteorological  tower.  This  swale  demon- 
strated an  effect  on  the  momentum  gradient  data,  as 
discussed  later.  A  complete  discussion  of  the  site,  in- 
cluding biota  and  climate,  is  given  by  Zeller  et  al.  (1989). 

Figure  3  is  a  photograph  of  the  physical  layout  of 
the  site,  which  includes:  a  10-meter  meteorology  tower 


Figure  3 — Pawnee  site  meteorological  tower  and  instrument 

shelter. 

equipped  with  micrometeorological  sensors  at  8  and  3 
meters;  an  instrument  shelter  located  21  m  southwest 
(21°  true  north  (TN)  azimuth)  of  the  tower;  fast-  and 
slow-response  chemical  ozone  sensors;  a  personal  com- 
puter; data  acquisition  system;  and  the  plumbing  and 
wiring  necessary  to  sample  the  air  at  the  tower  adjacent 
to  the  meteorological  sensors.  A  schematic  of  the  site  is 
given  by  Zeller  et  al.  (1989)  (fig.  1).  Table  2  lists  the 
major  instrumentation  used  to  collect  the  data  used  in 
this  analysis.  The  following  sections  present  further 
discussions  of  these  instruments,  and  possible  chemical 
and  temperature  effects  on  ozone  measurements.  Data 
acquisition  and  real-time  computations  made  during 
data  acquisition  are  discussed  in  appendix  A.  Measure 
spectra  and  cospectra  are  presented  in  appendix  B. 

The  Pawnee  site  was  operated  on  a  weekly  Tuesday 
morning  to  Tuesday  morning  basis  to  correspond  with 
other  acid  deposition  related  programs  and  experiments 


^5J65 


Figure  2 — Topographic  map  of  Pawnee  site  area  and  Pawnee  site  (*)  location. 
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Table  2.  —  Instrumentation  used  at  Pawnee  site. 


Instrument  or 
model 

Type 

Quantity 

Response 
time 

Range 

Denver  University 
CAAM  (1  &  2) 

Chemiluminescence 

Ozone  flux 

0.5  s 

0.2  to  400  ppb 

TECO  49  (1  &  2) 

UV  adsorption 

Ozone  concentration, 
ozone  gradient 

2  min 

1 .0  to  500  ppb 

AIR,  Inc. 

Model  FT-1AO-T 

Platinum  resistance 

Temperature 

0.01  s 
(4  m  S"1  wind) 

-50  to  +50°C 

Atmos.  Tech.  Inc. 
Model  BH-478/B3 

3-axis  sonic 

Wind  velocity 
momentum  flux 

0.05  s 

0  to  20  m  s1 

Compaq  386/20 

PC  computer 
(80386) 

Data  collection 
and  processing 

Data  Translations,  Inc. 

DT2805  board 

A/D  conversion 

Oto  10  v 

Thermocouple 

copper-constantan 

Temperature 
gradient 

±  10°C 

involving  weekly  averaged  samples  requiring  adher- 
ence to  national  acid-deposition,  data-collection  proto- 
cols. Therefore,  weekly  calibrations,  routine  mainte- 
nance, sensor  chemical  reagent  changes,  and  data  archi- 
val were  always  accomplished  on  Tuesday  mornings. 


Meteorological  Measurements 

The  eddy-correlation  approach  to  measuring  the 
vertical  fluxes  of  momentum,  heat,  and  mass  requires 
simultaneous  measurements  of  vertical  velocity  fluc- 
tuations, sensible  heat  fluctuations,  and  chemical  mass 
fluctuations.  Figure  3  shows  the  eddy-correlation  in- 
strumentation at  8  and  3  meters.  The  fast-response 
measurements  required  to  capture  the  eddy  fluctuations 
were  routinely  made  at  14  ±  1  Hz.  Wind  sensors,  tem- 
perature sensors,  and  chemical  sensors  must  all  have 
short  response  times  to  respond  at  this  sampling  rate, 
which  requires  very  sensitive  transducers.  All  sensors, 
however,  have  measurement  limitations  that  must  be 
considered  and  accounted  for.  This  is  particularly  true 
for  eddy-correlation  instruments  (Businger  1986)  be- 
cause the  ability  to  measure  eddy  fluctuations  is  limited 
by  sensor  response.  Specific  corrections  for  measured 
limitations  that  can  be  applied  to  eddy-correlation  flux 
data  are  presented  in  appendix  A  and  elsewhere  (Zeller 
et  al.  1989). 

Scalar  gradient  data  were  acquired  using  slower 
response  sensors.  Temperature-gradient  measurements 
were  made  with  3-mil  copper-constantan  thermocouples 
connected  in  series.  Solar  heating  of  the  thermocouples 
was  accounted  for  in  the  results.  Ozone  gradient  mea- 
surement are  discussed  in  appendix  A.  Wind  gradients 
were  obtained  from  average  sonic-anemometer  data. 

It  is  appropriate  to  investigate  the  response  of  the 
meteorological  and  chemical  sensors  as  a  function  of 


frequency  in  order  to  evaluate  the  ability  of  each  sensor/ 
data-logging  combination  to  detect  the  full  range  of 
energy-containing  eddies.  This  is  accomplished  by  trans- 
forming a  time  series  data  record,  e.g.,  the  vertical 
velocity  component  w(t)  or  any  other  measure  q(t)  as  a 
function  of  real  time  to  it's  power  spectrum,  a  function 
of  frequency.  In  this  way  the  measured  variance  density 
at  each  frequency  can  be  examined  and  compared  to  M- 
0  similarity  predictions: 


q'2=  |Sq(/)d/=|/Sq(/)dln/ 

where 


oo  oo 


[19] 


S{fl 

f 

q 


special  density; 
frequency  (s"1);  and 
measured  quantity. 


Equation  [19]  assumes  that  a  temporal  record  can  be 
used  to  represent  the  spatial  "size"  of  turbulent  eddys. 
This  assumption,  that  means  and  variances  measured  in 
time  are  equal  to  those  measured  in  space,  is  Taylor's 
"frozen  wave"  hypothesis  and  is  valid  for  stationary- 
homogeneous  turbulence  (Panofsky  and  Dutton  1984). 

During  the  measurement  program,  half-hour  time 
series  records  of  instantaneous  or  "raw"  time  series  data 
were  collected  at  either  ca.  40  Hz  or  ca.  14  Hz. 
These  data  were  translated  to  the  frequency  domain 
using  the  fast  Fourier  transform  (FFT)  and  power-spec- 
trum Fortran  programs  provided  in  "Numerical  Reci- 
pes" (Press  et  al.  1986).  The  data-handling  approach, 
normalizations,  and  data  presentations  are  similar  to 
those  used  by  Kaimal  et  al.  (1972):  block  averaging  was 
applied  to  the  spectra  and  cospectra.  Individual  spectra, 
S{f],  are  typically  normalized  and  plotted  as/S(/)/q*  vs. 
fzl  u  on  a  log-log  plot,  where /z/u  isthenondimensional 
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frequency.  Monin-Obukhov  similarity  predicts  a  -5/3 
slope  in  the  inertial  subrange  for  both  velocity  and 
scalar  spectra  when  plotted  against  frequency.  This 
translates  to  a  -2/3  slope  for  fS{f)  plotted  against  fre- 
quency. Temperature,  vertical  velocity,  and  ozone  power 
spectra  are  presented  in  appendix  B  along  with  vertical 
velocity-temperature  cospectra  and  vertical  velocity- 
ozone  cospectra. 


Ozone  Measurements 

By  combining  [7]  and  [8b],  the  eddy  diffusivity  for 
ozone  can  be  determined  by  measuring  both  the  flux  by 
eddy  correlation  and  the  concurrent  ozone  gradient: 


-  w'c'  Az 


dz 


Ac 


[20] 


The  measurement  of  w'c'  by  eddy  correlation  requires  a 
fast-response  chemical  sensor.  The  Ac  =  c(8  m)  -  c(3  m) 
measurement  is  very  sensitive  to  both  instrument  accu- 
racy and  precision.  Ozone  measurements  at  the  Pawnee 
site  were  made  with  four  ozone  instruments  (see  table 
2),  two  fast-response  chemiluminescent  ambient  air 
monitors  (CAAM) ,  and  two  slow-response  Thermo  Elec- 
tron Corp.  (TECO)  model  49  UV  absorption  ambient  air 
monitors. 

Ozone  Analyzers 

The  TECO  49  ozone  analyzer  has  been  used  for  years 
in  environmental  monitoring  programs  and  is  known 
for  its  stability  (TECO  1986).  It  can  provide  a  very 
precise  and  relatively  accurate  continuous  ozone  record 
month  after  month.  It  has  been  accepted  by  the  U.S. 
Environmental  Protection  Agency  as  an  "equivalent 
method"  (EQOA-0880-047).  The  TECO  49  measures  the 
extinction  of  254  nm  wavelength  light  by  ozone.  Beer's 
law  with  corrections  for  pressure  and  temperature  is 
used  to  calculate  the  ozone  concentration  in  parts  per 
billion  (ppb).  The  TECO  measures  mixing  ratio  as  op- 
posed to  mass  per  volume. 

The  CAAM  ozone  analyzer  is  a  prototype  instru- 
ment developed  and  constructed  by  Dr.  Donald  Stedman 
of  the  University  of  Denver.  Its  response  is  linear  to  the 
number  of  ozone  molecules  in  the  sample  air  stream 
(i.e.,  mass  per  volume).  Due  to  its  sensitivity,  it  is  not  as 
accurate  over  long  time  periods  as  the  TECO  49;  how- 
ever, its  linear  response  gives  it  a  relative  precision  that 
allows  for  its  use  in  measuring  ozone  fluctuations.  The 
CAAM  uses  a  blue-sensitive  photomultiplier  tube  to 
detect  the  light  produced  by  the  reaction  of  ozone  with 
the  chemiluminescent  dye  eosin-y  mixed  in  solution 
with  water,  ethylene  glycol,  and  copper  nitrate.  For 


more  detail  on  the  CAAM,  see  Ray  et  al.  (1986),  and 
appendix  A  in  Zeller  et  al.  (1989). 

CAAM  Response 

The  square-wave  response  time  to  an  ambient-air  to 
clean-air  step  change  for  the  CAAM  analyzer,  assuming 
a  first-order  sensor,  was  measured  in  the  laboratory  to  be 
about  0.14  second.  In  situ  this  response  time  deterio- 
rated to  either  0.7  second  if  the  step  change  of  clean  air 
was  introduced  at  the  tower  end  of  the  20  m  intake  line, 
or  0.3  second  if  the  step  change  of  clean  air  was  intro- 
duced at  the  intake  manifold  just  prior  to  entering  the 
CAAM.  Figure  4  shows  both  these  response  times. 
Diagrams  and  details  are  provided  in  Zeller  et  al. 
(1989).  For  the  ozone  flux  correction  considerations,  a 
time  response  of  0.5  second  was  used  as  representative 
of  the  actual  response  time  seen  by  the  CAAM  as  the 
intake  line  is  actually  in  equilibrium  with  the  ambient 
air  during  routine  sampling  but  tends  to  adsorb  and 
desorb  during  zero  air  tests. 

CAAM  Lag 

The  sample  transfer  time  for  a  sample  of  air  to  travel 
from  the  tower  inlet  to  the  CAAM  was  ca.  2  seconds. 
This  transfer  time  depended  on  the  flow  rate  of  the  high- 
volume  pump  and  the  time-in-use  of  the  pump  filters. 
The  flow  rate  averaged  ca.  155  liters  per  minute,  which 
included  the  additional  time  from  the  intake  manifold 
inside  the  instrument  shelter  to  the  CAAM  instrument 
itself.  Since  the  measure  w'c'(t)  requires  simultaneous 
values  at  any  given  t  for  both  w'(t)  and  c'(t),  it  was 
important  to  know  the  lag  time  required  for  the  sampled 
air  to  be  detected.  The  lag  time  was  checked  for  both 
intake  lines  routinely  on  Tuesday  mornings  and  incor- 
porated into  the  flux  calculations.  The  lag  time  for  both 
instruments  remained  constant  within  ±0.05  seconds 
throughout  the  study.  The  initial  2  seconds  shown  on 
figure  4  is  the  lag  time,  as  the  data  recording  for  the  data 
plotted  commenced  with  the  injection  of  zero  air  at  the 
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Figure  A — CAAM  signal  response  to  square-wave  ozone  input- 
April  23, 1989  (13:48  MST) 
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tower.  The  dampening  of  turbulent  fluctuations  of 
ozone  in  the  intake  tube  was  found  to  be  insignificant 
relative  to  measurement  of  atmospheric  turbulent 
fluxes  (Massman  1991). 

Ozone  Calibrations 

The  TECO  49  ozone  analyzers  were  calibrated  on  a 
weekly  basis  with  a  primary  standard  ozone  analyzer 
maintained  by  the  USDA  Forest  Service,  Rocky  Moun- 
tain Range  and  Experiment  Station,  strictly  for  the 
purpose  of  ozone  calibrations.  The  local  Forest  Service 
primary  standard  instrument  is  a  TECO  model  49PS. 
This  instrument  is  traceable  to,  and  compared  yearly, 
with  an  official  EPA  ozone  standard. 

The  CAAM  ozone  analyzers  were  calibrated  con- 
tinuously against  the  TECO  analyzers.  For  each  half- 
hour  period,  the  average  CAAM  output  signal  was 
compared  with  the  average  output  of  the  TECO  sam- 
pling the  same  air  arriving  through  the  same  intake 
manifold.  Figure  5  is  an  example  of  the  data  scatter 
between  the  half-hour  average  TECOl  and  CAAM2 
outputs  for  the  1-week  period  between  Tuesday  morn- 
ing June  27,  1989  and  Tuesday  morning  July  4,  1989. 
The  two  TECO  49  and  two  CAAM  ozone  analyzers  were 
named  TECOl,  TEC02,  CAAMl,  and  CAAM2  at  the 
Pawnee  site  for  calibration  and  maintenance  purposes. 
The  CAAM  instruments  were  calibrated  to  the  TECO 
instruments  on  a  half-hour  to  half-hour  basis  because 
the  TECOs  were  very  stable  over  a  1-week  period 
whereas  the  CAAMs  were  not  as  stable  over  the  same 
period  (note  there  is  some  scatter  in  fig.  5).  The  CAAM/ 
TECO  calibrations  were  done  after  TECOl  and  TEC02 
analyzers  were  aligned  with  each  other  to  provide  for 
accurate  dc/dz  ozone  gradient  data.  The  alignment 
procedure  is  described  in  appendix  A, 

Figure  6  is  an  example  of  a  nine-point  calibration  of 
a  CAAM  with  a  TECO  when  each  instrument  was  com- 
pared against  known  ozone  concentrations  under  con- 
trolled conditions.  The  figure  demonstrates  the  linear 
response  of  the  CAAM  instrument;  however,  the  cali- 
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TEC01  output  (ppb) 

Figure  5— CAAM2  (volt)  averages  vs.  TEC01  (ppb)  averages  for 
the  week  June  27-July  4, 1989. 


o 
S 

Q. 
3 
O 
CM 


CAAM2  =  (0.102  ±  0.001)TECO  2  +  (0  051  +  0.120) 
rho  =  0.9989 


10 


30  50  70 

TEC02  output  (ppb) 


90 


110 


Figure  6 — CAAM  vs.  TECO  9-point  calibration— April  5, 1989 

(12:30  MST). 

bration  shown  would  not  necessarily  hold  over  a  1-week 
period  even  though  the  CAAM  response  would  remain 
linear  at  all  times.  When  challenged  with  a  constant 
ozone  calibration  gas,  the  CAAM  instrument  outputs  a 
white-noise  signal  that  varies  ±1.7%  of  the  expected 
value.  For  instance,  for  a  30-ppb  ozone  calibration  gas 
the  CAAM  output  range  would  fluctuate  between  29.5 
and  30.5  ppb. 

TECO  Response  vs.  Ambient  Temperature 

One  of  the  more  interesting  results  of  this  study  is 
the  similarity  in  variation  with  time  of  the  Ac(t)  and 
A6[t)  measurements  between  8  and  3  meters.  Since 
A  c(t)  is  measured  by  subtracting  the  averages  of  the  two 
TECO  analyzers,  one  might  suspect  that  the  TECOs  are 
responding  to  temperature  as  well  as  to  ozone.  Although 
the  TECO  49  has  an  internal  adjustment  for  inlet  air 
temperature  (TECO  1986)  it  is  appropriate  to  verify  its 
response  to  varying  inlet  temperatures.  Figure  7  is 
a  plot  of  inlet  concentration  divided  by  the  TECO  output 
concentration  while  the  TECO  is  challenged  with  a 
continuous  stream  of  20-ppb  ozone  calibration  gas  inlet 
air  during  the  process  of  heating  the  inlet  air  from  25°C 
to  51°  C  over  a  period  of  20  minutes.  Additional 
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Figure  7— TECO  response  to  varying  inlet  air  temperature  over  20- 
minute  period— September  27, 1989. 
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temperature  response  tests  were  made  by  challenging 
the  TECO  with  80-ppb  inlet  air  at  22°  C  continuously  for 
10  minutes  and  again  at  51°  C  continuously  for  10 
minutes  with  the  following  results: 

Inlet  Calibration  gas  TECO  49 

temperature  (°  C)     input  c  (ppb  03)      output  c  (ppb  03) 


22.1 
51.4 


80.0 
79.9 


79.8 
80.8 


It  was  concluded  that  the  TECO  49  analyzer  is  not 
sensitive  to  ambient  temperature. 

Ozone  Depletion  Due  to  NOx 

Ozone  is  a  very  reactive  gas.  The  assumption  that  no 
ozone  is  lost  to  chemical  reactions  within  the  layer  of  air 
between  3  and  8  m  is  questionable  (Lenschow  1982). 
Ozone  is  a  member  of  the  03-NO-N02  triad  that  takes 
part  in  the  following-reactions: 


N02  +  hv 
O  +  02- O 
NO  +  O. 


NO  +  O 


3 

NO 


[21] 


2  +  o2 


The  first  and  third  reactions  are  much  slower  than  the 
second  and  are  therefore  of  primary  concern  (Lenschow 
and  Delany  1987). 

Models  that  calculate  ozone  flux  based  on  reactions 
with  NOx  [21]  and  without  NOx  give  similar  results  and 
compare  well  with  measured  data  if  certain  criteria  are 
met  (Fitzjarrald  and  Lenschow  1983).  Fitzjarrald  and 
Leschow  (1983)  give  the  criteria 


NO 


<  0.1 


[22] 


NOx 


(where  NOx  is  NO  +  N02  +  other  oxides  of  nitrogen)  as 
a  limit  to  the  importance  of  incorporating  [21]  into 
ozone  flux  calculations  below  2  m.  In  effect,  [22]  is  a 
criterion  for  disregarding  a  net  change  (production  or 
loss)  of  ozone  due  to  [21].  Based  on  the  fact  that  ozone 
concentrations  are  typically  an  order  of  magnitude  greater 
than  NOx  concentration,  Lenschow  and  Delany  (1987) 
state  that  any  error  in  triad  03-NO-N02  flux  measure- 
ments is  least  serious  for  ozone,  particularly  in  a  non- 
urban  environment.  They  claim  flux  errors  should  be 
<10%.  Concentrations  of  NO  and  NOx  were  measured, 
albeit  not  simultaneously,  at  the  Pawnee  site  during  the 
summer  of  1988  (Stocker  et  al.  1989).  The  average  NO 
concentration  during  July  19  through  August  2,  1988 
was  0.6  ppb  ±  0.4  ppb.  The  average  NOx  concentration 
during  June  21-28,  1988  was  5  ppb  ±  3  ppb.  Stocker  et 
al.  (1989)  also  noted  that  higher  values  of  NO  and  N02 
flux  were  seen  up  to  8  hours  following  rainfall  events. 
Concentrations  for  both  pollutants  vary  diurnally  and 


were  not  measured  simultaneously  as  the  same  instru- 
ment is  used  to  measure  both  NOx  and  NO  (Zeller  1988; 
appendix  A).  Using  the  above-measured  average  con- 
centrations for  characteristic  values  of  NO  and  NOx  at 
the  Pawnee  site  provides  an  order-of-magnitude  ratio 
[22]  of:  0.6/5  =  0.12  (12%)  for  NO  compared  to  NOx. 
Therefore,  the  chemical  reaction  of  ozone  with  NO  is  a 
possible  source  of  error;  however,  these  measurements 
were  taken  above  2  m  where  it  is  most  likely  that  an 
equilibrium  in  the  03-NO-N02  triad  exists  (Lenschow 
and  Delany  1987).  In  addition  to  NO,  aerosols  sus- 
pended in  the  atmosphere  will  also  reduce  ozone  con- 
centrations. 

Nonmethane  hydrocarbons  (NMHC)  also  play  a  role 
in  the  production  of  ozone,  particularly  for  NOx  con- 
centrations above  1.5  ppb  (Liu  et  al.1987)  which  have 
been  measured  at  the  Pawnee  site.  The  emission  of 
NMHC  at  the  Pawnee  site  is  unknown;  however,  it  is 
expected  to  be  small  (Williams  1990)  and  not  a  large 
contributor  to  ozone  production. 

For  the  purpose  of  this  paper,  it  is  assumed  that 
ozone  approximates  a  nonreactive  trace  gas,  although 
this  assumption  is  actually  dependent  upon  NOx  con- 
centrations. 


RESULTS  AND  DISCUSSION 

Based  on  completeness  and  continuity  of  the  quan- 
tities measured  and  the  operating  integrity  of  the  overall 
data-collection  system,  the  data  periods  edited  and 
analyzed  are  May  16-23, 1989  (Julian  dates  136-143)  and 
June  6  through  July  11,  1989  (Julian  dates  157-192). 
Appendix  C  shows  the  crucial  data  records — fluxes, 
averages,  and  gradients  of  ozone,  temperature,  and  wind 
speed  along  with  net  and  total  radiation — plotted  in 
their  entirety  (unedited).  The  diurnal  nature  of  these 
measures  is  evident.  The  data  in  appendix  C  represent 
approximately  1,900  half-hour  sampling  periods.  After 
applying  the  editing  criteria  presented  in  appendix  A, 
table  A-2,  approximately  730  data  records  were  avail- 
able to  provide  for  the  micrometeorological  analyses 
presented  below.  All  the  data  presented  below  are 
edited  and  all  the  ozone  flux  data  are  corrected  for  heat 
flux,  as  discussed  in  appendix  A.  Other  flux  corrections 
discussed  in  appendix  A  and  in  Zeller  et  al.  (1989)  are 
tested  and  discussed  below,  but  they  are  not  included  in 
the  determinations  of  (J>q's  nor  Kq's.  Reference  to  edited, 
uncorrected  data  will,  therefore,  mean  that  they  were 
corrected  for  heat  flux  as  just  discussed;  however,  they 
were  not  corrected  for  frequency  response  or  the  other 
response-related  flux  corrections  just  alluded  to. 

Several  of  the  following  discussions  present  data 
distributions  as  functions  of  £  =  z/L.  The  value  of  £  was 
calculated  from  the  data  using  L  from  [13],  with  the 
values  of  u*  and  w'0'  provided  by  averaging  the  mea- 
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sured  values  at  3  and  8  meters.2  Because  the  flux  data 
from  two  heights  are  averaged,  an  applicable  average 
height  must  be  chosen.  In  as  much  as  wind  distributions 
in  the  surface  layer  are  theoretically  logarithmic  with 
height  during  neutral  conditions  (Businger  1986),  the 
elevation  for  which  calculated  t,  and  K„  values 
presented  herein  are  applicable  is  4.9  meters,  the  geo- 
metric mean  measurement  height:  (3  x  8)1/2. 

Gradients  of  Ozone,  Temperature  and  Horizontal 
Wind  Speed 

The  approach  used  to  obtain  gradients  of  ozone, 
temperature,  and  wind  speed  is  discussed  in  appendix 
A.  The  structures  of  the  simultaneous  time-series  for 
ozone-  and  temperature-gradient  measurements  are 
similar  and  appear  to  track  each  other,  as  shown  in  the 
appendix  C  plots  of  the  ozone  and  temperature  gradi- 
ents and  also  in  Zeller  et  al.  (1990,  fig.  1).  Care  was  taken 
to  ensure  that  the  TECO  sensors  were  monitoring  ozone 
and  they  were  not  sensitive  to  temperature  (see  appen- 
dix A).  To  demonstrate  the  variation  with  time  of  each 
measure,  figure  8  shows  the  gradients  of  ozone  and 
temperature  for  the  period  June  13-21,1989  (Julian  dates 
164-172).  Note  that  although  the  two  gradients  track  in 
general,  they  do  not  track  exactly.  Both  ozone  and 
temperature  gradients  are  sensitive  to  cloud  cover 
which  can  be  inferred  by  examining  the  solar-radiation 
data  corresponding  to  the  periods  when  the  extent  of  the 
day/night  gradient  cycles  are  minimal  (e.g.,  Julian  days 
164  and  168  compared  to  days  165, 166, 167,  and  169  in 
fig.  8).  The  Julian  day  plotted  in  all  figures  represents 
midnight  beginning  the  new  Julian  day;  for  example,  the 
data  for  June  17,  1989  starts  at  Julian  day  168  and  ends 
at  Julian  day  169  on  figure  8.  This  convention  is  fol- 
lowed in  all  figures  and  in  appendix  C. 

The  ozone-gradient  data  appear  to  reach  a  minimum 

each  morning  at  7:30  MST  (Mountain  Standard  Time) 

and  increase  slowly  with  time  of  day  until  early  evening, 

17:30  MST,  when  they  increase  rapidly.  Figure  9,  a 

diurnal  plot  of  the  mean,  and  variance  ozone-gradient 

data  for  the  entire  period,  shows  this  behavior  in  detail. 

On  average,  the  ozone-gradient  value  falls  rapidly  from 

a  nighttime  value  of  ca.  0.6  ppb  m"1  to  0.16  ppb  m"1 

between  sunrise  (-4:30  MST)  and  7:30  MST,  then  slowly 

increases  to  ca.  0.2  ppb  m"1  by  16:00  MST,  levels  off,  and 

increases  rapidly  after  17:30  MST  as  sunset  approaches. 

Ozone  concentrations  in  the  surface  layer  behave  in  the 

reverse  (fig.  10):  the  concentration  of  ozone  is  lowest  at 

night,  increases  at  sunrise  to  a  maximum  at  midday, 

then  decreases  with  sunset.  Ozone  is  known  to  pool 

above  the  surface  layer  at  night  when  vertical  transport 

is  minimal  (ozone  near  the  surface  is  lost  to  deposition 

2Both  Kaimal  (1989)  and  Wyngaard  (1989a)  agreed  that  this  averaging 
approach  is  reasonable  for  applying  the  information  available 
from  both  levels. 
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Figure  8 — Ozone  and  temperature  gradient  data  for  June  1 3-21 , 1 989. 
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Figure  9 — Mean  and  variance  diurnal  ozone-gradient  data  for  the 
entire  period:  May  16-23  and  June  7-11, 1989. 


100 


90  - 


«-  80 
2  70 


Q. 

E 

a> 


n 

0) 

c 
o 

N 

O 


60  - 


50 


40 


30 


20 


c  8  meters 
c  3  meters 
0  8  meters 
0  3  meters 


10 


164 


166 


168 

Julian  date 


170 


172 


Figure  10 — Ozone  concentration  and  temperature  at  3  and  8 
meters  for  June  13-20, 1989. 

and  chemical  reactions  with  surface  elements),  then 
mix  rapidly  downward  during  daylight  with  the  eradi- 
cation of  the  nighttime  inversion  as  a  result  of  surface 
solar  heating  (Edinger  1973;  Vukovich  etal.  1985;  Zeller 
et  al.  1977).  This  phenomenon  explains  the  average 
values  presented  in  figure  9:  initial  morning  mixing  of 
air  from  the  surface  up  is  rapid  and  the  source  of  ozone 
is  limited  vertically,  hence  the  vertical  gradient  is  rap- 
idly reduced  to  a  minimum.  As  the  inversion  base 
continues  to  lift,  additional  ozone  from  aloft  or  possibly 
generated  in  situ  is  available  for  downward  transport. 
The  gradient  increases  until  an  equilibrium  is  reached 
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in  the  late  afternoon  between  regional  ozone  produc- 
tion, advection,  and  surface  ozone  deposition. 

Figure  11  plots  ozone  gradient  data  vs.  temperature 
gradient  data.  Negative  temperature  gradients  are  asso- 
ciated with  daylight  hours  and  positive  temperature 
gradients  with  nighttime  hours.  Note  that  in  general 
there  is  no  correlation  between  the  two  gradients  during 
the  day,  but  at  night  there  is  some  correlation.  The 
nighttime  correlation  line  plotted  is  A  03  =  (1.62  ±  0.13) 
A  T  +  (0.17  ±  0.13)  (±  values  are  standard  errors  for  slope 
and  intercept,  respectively)  with  a  correlation  coeffi- 
cient (rho)  of  0.65.  Strong  positive  temperature  gradi- 
ents are  indicative  of  limited  mixing,  which  accounts 
for  limited  ozone  turbulent  transport  and  hence  larger 
ozone  gradients. 


Variance  of  Ozone  Temperature,  and  Vertical 
Wind  Speed 

The  variances  of  scalars  and  vertical  velocity  appear 
to  obey  M-0  scaling  in  the  unstable  surface  layer  better 
than  fluxes  and  gradients  (Panofsky  and  Dutton  1984). 
With  the  goal  of  investigating  the  similarity  of  Kq's,  it  is 
worthwhile  to  verify  the  behavior  of  the  Pawnee  site 
data  variances  for  w,  6,  and  c  with  the  empirical  rela- 
tionships provided  by  other  researchers.  Equations  [23a] 
and  [23b]  are  two  empirical  relationships  for  vertical 
velocity  and  scalar  quantities,  respectively  (Panofsky 
and  Dutton  1984): 


w 


1.25(1  -3Q1/3     for  ^<0 


u 


q    =   5(l-16g_1/2     for  £<0 


[23a] 


[23b] 


where    cQ    =  (q,2)1/2 


scalar    quantities,  and 


/2  a/2 


1.0 


(Day) 


(Night) 


Ac  =  (1.62  ±  0.13)A5  +  (0.17  ±  0.13) 
rho  =  0.65 


-0.1  0.1  0.3  0.5 

Temperature  gradient  (°C  m-1) 

Figure  11 — Ozone  gradient  vs.  temperature  gradient  for  the  entire 
data  collection  period. 
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Figure  12 — Nondimensional  ozone  variance  <jc  I  Ic.l  vs.  stability 
parameter  (z/L)  for  (a)  3-meter  and  (b)  8-meter  data — May  16-July 
11, 1989. 

The  constants  in  these  equations  are  not  particu- 
larly well  established  and  vary  from  investigator  to 
investigator.  The  1/3  and  -1/2  powers,  however,  do 
describe  typical  slopes  of  measured  variances  versus  £ 
for  £  <  -0.5,  the  local  free-convection  regime  (Wyngaard 


1989b).  Figure  12  shows  oc  / 


for  ozone  variance  for 


3  and  8  meters  as  measured  by  the  CAAM  ozone  sensors; 
the  theoretical  plots  of  [23b]  are  adjusted  for  z  =  3  and 
z  =  8  meters,  respectively.  The  data  scatter,  less  at  8 
meters  than  at  3  meters,  is  typical  of  that  experienced  by 
others  (Panofsky  and  Dutton  1984);  also,  the  trend  with 
unstable  £  follows  the  -1/2  slope  (solid  line)  closely. 

Figure  13  shows  a  el  ^  |  for  temperature.  Again  the 
data  scatter  at  3  meters  is  greater  than  that  at  8  meters 
and  the  data  slope,  -1/2,  follows  the  theory  closely 
(Panofsky  and  Dutton  1984).  The  relative  vertical  loca- 
tion of  the  theoretical  line  is  a  function  of  the  constants 
used  in  [23b].  This  investigation  does  not  attempt  to 
establish  constants  specific  to  the  Pawnee  site  for  [23a] 
and  [23b]. 

Figure  14  shows  o\.r  /  u  for  momentum.  For  this 
quantity  the  measured  data  scatter  at  8  meters  is  greater 
than  at  3  meters.  Vertical  turbulence  is  a  function  of  z; 
hence  the  spectrum  of  eddies  contributing  to  aw  at  8 
meters  is  broader  than  the  spectrum  of  eddys  contribut- 
ing to  o~w  at  3  meters.  The  data  follow  the  1/3  slope  of 
the  theoretical  curves,  [23a],  as  plotted. 
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Figure  13 — Nondimensional  temperature  variance  v9l  \t).\  vs.  stabil- 
ity parameter  (z/L)  for  (a)  3-meter  and  (b)  8-meter  data — May  1 6- 
July  11, 1989. 


Figures  12-14  demonstrate  that  individual  eddy 
measurements  for  ozone,  temperature  and  momentum 
at  3  and  8  meters  are  reasonable  and  follow  M-0  scaling 
within  the  bounds  typically  experienced  by  microme- 
teorological  researchers. 


Surface  Roughness 

The  determination  of  velocity,  temperature,  and 
mass  profiles  in  the  surface  layer  requires  a  knowledge 
of  the  surface  roughness,  z0,  and  the  zero-plane  dis- 
placement. For  the  Pawnee  site,  the  zero-plane  dis- 
placement is  essentially  zero  due  to  the  shortgrass 
prairie  canopy  environment  (individual  grass  element 
heights  ca.  25  cm);  hence  [17a-c]  do  not  include  it.  It  is 
appropriate  to  limit  the  determination  of  z0  to  near- 
neutral-stability  data  where  l£l  <  0.03  (eg.,  Fazu  and 
Schwerdtfeger  1989).  Figure  15  shows  the  results  of 
133  zQ  calculations  meeting  this  condition  plotted 
against  their  corresponding  wind  speeds  at  8  meters. 
The  value  z0  was  calculated  by  combining  [17a]  with 
[18a]  and  [18b].  In  this  way  zQ  is  based  on  measured  u* 
not  A  u.  The  average  zQ  for  the  Pawnee  site,  for  the  data 
period  analyzed,  is  1.60  cm  with  standard  deviation  of 
0.89  cm.  This  value  falls  within  the  expected  value  for 
"fairly  level  grass  plains"  (Arya  1988). 
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Figure  14 — Nondimensional  vertical-velocity  variance  <rw  /  u,  vs. 
stability  parameter  (z/L)  for  (a)  3-meter  and  (b)  8-meter  data — May 
16-July  11, 1989. 


Fluxes  of  Ozone,  Sensible  Heat,  and  Momentum 


The  unedited  and  uncorrected  fluxes  of  ozone,  w'c' 
(ppbms"1)  sensible  heat,  w'r3'(°Cms"1)  and  momentum, 
u'w'  (m2  s"2)  are  presented  in  appendix  C.  This  section 
presents  edited-corrected  and  edited-uncorrected  fluxes 
measured  at  3  and  8  meters. 

The  concept  of  a  constant-flux  layer  is  often  con- 
fused with  the  idea  that  the  flux  of  a  given  quantity  at 
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Figure  1 5 — Surface  roughness  zo  vs.  8-meter  wind  speed  for  K '<0.03. 
Average  zo=  1 .60  cm  standard  deviation  =  0.89  cm.  — May  1 6-July 
11, 1989. 
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two  separate  elevations  must  be  the  same  (Wyngaard 
1989a).  Rather  it  is  the  flux  divergence,  3w'q'  /  dz ,  that 
should  be  constant  with  height  in  the  surface  and 
convective  boundary  layers  (Wyngaard  1989b).  This 
implies,  for  instance,  that  the  flux,  w'q',  would  be  a 
maximum  for  heat  and  momentum  at  the  earth's  surface 
or  source  of  these  fluxes,  and  decrease  linearly  to  ap- 
proximately zero  at  the  top  of  the  boundary  layer.  On  the 
other  hand,  the  source  of  ozone  is  above  the  surface 
layer  (at  least  initially  in  the  morning)  in  non-ozone- 
producing  rural  locations  like  the  Pawnee  site.  For 
ozone  the  flux,  w'c',  would  be  a  maximum  at  some 
elevated  location  experiencing  maximum  ozone  con- 
centrations (i.e.,  its  source  for  vertical  transport),  and 
decrease  linearly  to  the  surface  where  it  doesn't  quite 
become  zero  but  is  continually  deposited  at  the  surface. 
One  would,  therefore,  expect  measured  u'w'  and  w'0' 
to  be  greater  at  3  meters  compared  to  8  (Wyngaard 
1989b),  and  measured  w'c'  for  ozone  to  be  greater  at  8 
meters  compared  to  3  based  on  the  above  argument. 

Haugen  et  al.  (1971)  have  used  the  criterion  ±  20% 
to  declare  a  given  flux  constant  within  the  surface  layer. 
In  practice,  daytime  boundary  layers  exceed  1  km; 
therefore,  the  distance  between  8  and  3  meters  relative 
to  the  top  of  the  boundary  layer  would  not  be  enough  to 
detect  flux  differences  given  the  limits  of  measurement 
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Figure  1 6 — Comparisons  of  the  uncorrected  fluxes  of  (a)  momentum 
(Reynolds  stress)  and  (b)  sensible  heat,  as  measured  at  8  and  3 
meters— May  1 6-July  11,1 989. 


accuracy  and  all  flux  measurements  between  3  and  8 
meters  ought  to  be  equal  statistically. 

Figure  16  compares  the  uncorrected  fluxes  of  (a) 
momentum  (Reynolds  stress)  and  (b)  sensible  heat  as 
measured  at  8  and  3  meters.  Although  in  both  figures  the 
measured  fluxes  at  3  meters  are  less  than  the  measured 
fluxes  at  8  meters  contrary  to  the  above  argument,  the 
statistical  differences  are  not  significant  as  the  95% 
confidence  intervals  fall  within  the  range  where  they 
are  equal.  For  a  subset  of  this  uncorrected  momentum- 
flux  data,  the  3-meter  flux  tended  to  be  slightly  greater 
(see  discussion  below  in  this  section). 

The  differences  between  the  8-  and  3-meter 
uncorrected  fluxes  of  ozone,  as  shown  in  figure  17  do 
not  behave  the  same  as  heat  and  momentum.  The  range 
for  constant  ozone  flux  at  both  heights  falls  at  the  edge, 
practically  outside,  of  the  95%  confidence  interval, 
which  is  statistically  significant.  That  is,  there  is  95% 
statistical  confidence  that  the  1:1  ratio  line  falls  outside 
the  measured  data.  Based  on  the  constant  flux  diver- 
gence argument,  the  flux  of  ozone  in  particular  should 
be  a  little  greater  at  8  compared  to  3  meters  but  it 
shouldn't  be  as  discernible  in  the  measurements  made. 
It  is  possible  that  ozone  is  depleted  chemically  at  3 
meters  faster  than  at  8  meters  because  of  surface  ozone 
sinks:  surface  canopy,  soil,  and  nitrogen  species.  How- 
ever, as  presented  previously,  ozone  should  be  in  chemi- 
cal equilibrium  above  2  meters  (Lenschow  and  Delany 
1987). 

One  possibility,  emissions  of  NO  (see  [21])  at  the 
Pawnee  site,  is  somewhat  dependent  on  soil  moisture. 
The  greater  the  soil  moisture  experienced  follow- 
ing rain  events,  the  more  NO  emitted  from  the  surface. 
Emissions  of  NO  are  neither  constant  with  time  nor 
diurnal  in  variability  (Stocker  et  al.  1989);  therefore, 
figure  17  should  show  even  more  scatter,  particularly 
about  the  1:1  ratio  if  NO  emissions  were  responsible, 
because  the  data  represent  both  moist  and  dry  periods. 
Another  possibility,  based  on  the  analysis  to  follow,  is 
that  the  known  errors  in  flux  measurements  as  dis- 
cussed in  appendix  A  and  Zeller  et  al.  (1989)  may 
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Figure  17— Ozone  flux  at  3  meters  vs.  ozone  flux  at  8  meters— May 

16-July  11, 1989. 
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account  for  the  measured  difference  in  ozone  flux  be- 
tween the  two  heights. 

Figure  18  compares  the  corrected  ozone  flux  at  the 
two  measurement  heights  for  two  weeks  of  data:  May  16- 
23  and  June  6-13,  1989,  subsets  of  the  total  data  base 
analyzed.  (Corrections  were  limited  to  this  2-week  test 
period  because  the  corrections  are  considerably  compli- 
cated to  make,  they  are  based  on  theoretical  cospectra, 
and  they  add  uncertainty  to  the  data.)  The  figure  dem- 
onstrates that  the  flux  corrections  may  account  for  the 
difference  in  measured  ozone  flux.  The  95%  confidence 
interval,  although  considerably  larger,  now  includes 
the  range  (1:1  ratio  line)  where  the  fluxes  would  be  equal 
statistically.  The  corrected  fluxes  have  been  plotted 
indicating  stability  because  corrections  for  stable  ozone 
fluxes  tend  to  be  very  large;  greater  than  10  times  the 
original  data  in  some  cases  (Zeller  et  al.  1989).  In 
general,  the  corrected  data  (fig.  18)  are  broadly  scattered 
compared  to  the  uncorrected  data  (fig.  17).  On  average, 
the  corrected  flux  at  8  meters  is  1.45  times  greater  than 
the  original  uncorrected  flux  data,  and  1 .84  times  greater 
than  the  original  uncorrected  data  for  3  meters. 

Figure  19  shows  the  uncorrected  residuals  for 
uncorrected  ozone  flux  [(w'c'  (8  m)  -  w'c'  (3  m)]  of 
figure  17  plotted  against  time  of  day.  The  data  are  too 
scattered  to  imply  much  except  that  the  differences  tend 
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Figure  18 — Corrected  ozone  flux  at  3  meters  vs.  corrected  ozone 
flux  at  8  meters— May  16-23  and  June  6-13, 1989. 


to  be  slightly  greater  during  midday  when  the  maximum 
ozone  flux  is  occurring.  This  tends  to  support  the  need 
for  flux  corrections,  particularly  with  slower  response 
sensors  sited  at  elevations  where  frequency  response  of 
the  sensors  are  inadequate.  The  larger  correction  values 
for  stable  cases  are  worrisome  however. 

Flux  corrections  were  also  made  on  the  momentum 
flux  and  sensible  heat  flux  data  for  the  same  data  subset: 
May  16-23  and  June  6-13,  1989  (fig.  20).  For  mo- 
mentum flux  (a)  the  corrections  are  not  as  severe  as 
those  for  ozone  flux.  Note  that  the  uncorrected  and 
corrected  momentum  flux  experienced  slightly  greater 
values  at  3  meters  than  at  8  meters  for  this  data  subset. 
The  95  %  confidence  interval  shown  in  figure  20a  is  only 
20%  larger  for  the  corrected  momentum  flux  compared 
to  uncorrected  flux,  considerably  less  than  the  300% 
increase  in  the  95%  confidence  interval  experienced  for 
ozone  flux  corrections.  On  the  average,  corrected  mo- 
mentum flux  for  both  heights  were  the  same:  1.34  times 
the  original  flux.  Sensible  heat-flux  corrections  (fig. 
20b)  behave  similarly  to  those  for  momentum  flux:  data 
scatter  as  reflected  in  the  data  and  in  the  95%  confi- 
dence interval  is  ca.  20%  greater;  corrections  for  8 
meters  are  1.21  times  the  uncorrected  flux  for  stable  and 
unstable  situations;  corrections  for  3  meters  are  1.21 
times  the  positive  (unstable)  flux  and  1.32  times  the 
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Figure  20a — Corrected  momentum  flux  at  3  meters  vs.  corrected 
momentum  flux  at  8  meters— May  16-23  and  June  6-13, 1989. 
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Figure  19 — Ozone  flux  residuals:  [(8-meter  flux)  -  (3-meter 
vs.  time  of  day — May  16-July  11, 1989. 
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sensible  heat  at  8  meters— May  16-23  and  June  6-13, 1989. 
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negative  (stable)  flux.  Momentum-flux  corrections  are 
slightly  greater  because  of  the  added  correction  for 
horizontal  wind  shadowing,  while  both  sensible-heat 
and  momentum-flux  corrections  are  less  than  ozone 
because  their  sensor  frequency  responses  are  consider- 
ably faster  than  those  of  the  ozone  CAAM  sensors. 


Dimensionless  Gradients  of  Ozone,  Sensible  Heat, 
and  Momentum 

As  discussed,  the  dimensionless  gradients  of  mo- 
mentum, sensible  heat,  and  mass,  [I5a-c],  are  expected 
to  be  universal  functions  of  £  according  to  M-O  similar- 
ity. Figure  21  shows  the  values  of  4>m,  4*^,  and  (}>c 
using  the  edited,  uncorrected  (see  first  paragraph  under 
"Results  and  Discussion")  Pawnee  site  data  plotted 
linearly  as  functions  on  £.  The  results  are  valid  for  z  =  4.9 
m  as  the  measured  data  from  both  3  and  8  m  are  used  to 
calculate  the  4>  values.  The  empirical  lines  annotated  in 
figure  21  are  [16a-b]  where  [166]  is  used  for  (}>c  in  figure 
21c  (i.e.,  <J)C  =  cj>h).  These  are  the  expressions  for  <J> 
established  by  Businger  (1971).  The  other  lines  are 
regression  models  of  the  data  using  £  as  the  independent 
variable  (discussed  below).  Table  3  presents  the  linear 
regression  results  obtained  by  comparing  the  measured 
values  of  4>:  cf>mm,  <t>km>  and  4>cm,  with  the  empirical  or 
theoretical  values  <}>mT,  <J>Jjt>  and  4>ct  (the  second  sub- 
script is  added  to  distinguish  between  measured,  m,  and 
theoretical,  T).  These  regression  results  are  included  to 
show  the  statistical  differences  between  the  Pawnee- 
site  measurement  and  those  established  for  smooth 
terrain.  The  regression  results  were  separated  by  wind 
direction  because  of  the  wind-direction  dependence  on 
<j>mm  shown  in  figure  21a.  Measured  values,  cj>mm,  are 
considerably  lower  for  wind  directions  with  azimuth 
250°  clockwise  through  azimuth  030°.  This  dramatic 
wind  direction  related  distribution  in  the  values  of  $mm 
is  surprising  as  the  z0  analysis  based  on  measured  u^ 
(fig.  15)  did  not  reveal  a  dependence  on  wind  direction. 
This  is  also  true  when  zQ  is  calculated  for  all  stabilities. 
There  are  two  possible  causes  for  the  dependence  of 
wind  direction  on  <f>m:  climatological — westerly 
through  northerly  wind  speeds  tend  to  be  greater  at  the 
Pawnee  site;  and/or  terrain  induced — as  described  pre- 
viously, there  is  a  drainage  swale  ca.  400  m  west  of  the 
site.  Measured  momentum  flux  tends  to  be  greater  for 
faster  wind  speeds  regardless  of  the  wind  direction, 
while  measured  wind  shear,  du/dz,  at  the  Pawnee  site 
tends  to  be  smaller  by  a  factor  of  two  for  westerly 
through  northerly  winds.  None  of  the_other  measured 
flux  nor  gradient  data,  (w'0',  w'c',  ddldz,  and  3c/3z) 
show  dependence  on  wind  direction.  The  fact  that  3  u/ 
dz  shows  a  dependence  on  wind  direction  and  that 
larger  u'w'  values  are  simply  related  to  faster  wind 
speeds  explains  why  the  terrain  effect  shows  up  for  4>m, 
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Figure  21a — Dimensionless  wind  shear  (momentum)  (<J>m)  data  re- 
sults— theoretical  dimensionless  wind  shear;  and  best  regres- 
sion model  wind  shear  for  east-  and  west-wind  categories  vs. 
stability  parameter  (z/L). 
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Figure  21  b — Dimensionless  sensible-heat  gradient  (4>h)  data  results; 
theoretical  sensible-heat  gradient;  and  regression-model  sen- 
sible-heat gradients  vs.  stability  parameter  (z/L). 
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Figure  21c — Dimensionless  ozone  gradient  (4>c)  data  results;  theo- 
retical sensible-heat  gradient;  and  regression-model  ozone  gra- 
dients vs.  stability  parameter  (z/L). 

[15],  and  not  for  z0,  [17].  The  data  results  specific  to 
momentum  will  hereafter  be  segregated  by  wind  direc- 
tion where  azimuth  250°  clockwise  through  030°  will  be 
referred  to  as  "west"  and  representative  of  "irregular  or 
rough  terrain,"  and  wind  direction  azimuth  030°  clock- 
wise through  250°  will  be  referred  to  as  "east"  and 
representative  of  "smooth  terrain." 

The  best  correlations  with  smooth  terrain  (table  3) 
for  <j>  values  are  for  unstable  <}>c,  unstable  <J>h  and  stable 
(}>m.  The  data  distribution  for  values  shown  in  figure  21 
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Linear  regression  4>qm  (measured)  vs.  4>qT 

(theory). 

%T 

No 

Corr. 

wind 
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btaDiliTy^ 
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points 
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0.72 
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0.69 
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0.37 
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0.71 
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2.06 
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0.38 
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0.36 
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s 

East 
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0.42 

s 

All 
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0.21 

0.63 

1.23 
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0.47 

s 

All 

2.16 

0.07 

0.00 

1.24 

231 

0.45 

u 

West 

3.35 

0.41 

-0.80 

0.66 

173 

0.53 

u 

Fact 

3.73 

0.20 

-0.76 

0.45 
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0.72 

u 

All 
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0.19 

-0.65 

0.54 
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0.63 

u 

All 
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0.05 

0.00 

0.57 
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0.58 

s 

West 

1.98 

0.30 

-0.39 

1.11 
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0.56 

s 

East 

1.08 

0.20 

0.30 

0.87 

114 

0.46 

s 

All 

1.32 

0.17 

0.15 

1.02 

214 

0.47 

s 

All 

1.43 

0.05 

0.00 

1.02 

214 

0.47 

vs  <t>mT 

u 
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0.77 

0.13 

-0.20 

0.19 

159 

0.43 

u 

West 

0.52 

0.02 

0.00 

0.19 

159 

0.41 

u 

East 

0.91 

0.17 

0.49 

0.36 

319 

0.29 

u 

East 

1.58 

0.03 

0.00 

0.37 
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0.19 

s 

West 

0.95 

0.13 

-0.44 

0.51 
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0.56 

s 
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0.64 

0.03 

0.00 

0.52 
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0.53 

s 

East 

1.13 

0.10 

0.67 

0.48 
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0.71 

s 

East 

1.53 

0.03 

0.00 

0.51 

131 

0.66 

+  u:  unstable;  s:  stable 

*  east:  030°  -  250°  azimuth  (clockwise) 

west:  250°-030°  azimuth  (clockwise) 

all:     all  wind  directions  combined. 
§4>cT=  <t>hT 

SE:     Standard  Error 

is  typical  in  the  literature  (e.g.,  Korrell  et  al.  1992;  Fazu 
and  Schwerdtfeger  1989)  where  reported  unstable  val- 
ues for  4>m  and  4>h  are  widely  scattered  with  values 
ranging  between  0.0  and  1.5.  Standard  errors  for  (J>qin 
comparisons  with  <J>^T  theory  as  given  in  table  3  are  not 
typically  presented  m  the  literature;  however,  the  rela- 
tively large  values  for  y-standard  errors  associated  with 
the  Pawnee  site  results  appear  to  be  normal  judging  by 
the  similar  ranges  in  published  data  scatter. 

Four  regression  models  with  £,  the  independent 
variable,  were  fitted  to  the  measured  4>  data  to  provide 
a  means  to  evaluate  eddy  diffusivities: 


qm 
qm 


=  ba  +  b2£ 
=  ba  +  b2  UYh- 


^qm  =  bl  +  eXP(b2  +  b3^ 
^qm  =  bl  +  exP(b 


2  +  b3^ 


t  >  o 

alH 


[24a] 
[24b] 
[24c] 
[24d] 


Here  bj  are  the  separate  parameter  constants  for  each 
model  and  for  each  comparison  (i.e.,  ba  [24a]  *b1  [24b]). 
Table  4  presents  the  regression  results  including  the 
standard  errors  on  each  parameter  and  the  correlation 
coefficient  for  each  model  result.  Regressions  were 
made  for  separated  east-  and  west-wind  categories  in 
addition  to  an  "all"  wind-direction  category.  The  corre- 


lation coefficients  are  similar  in  value  to  those  in  table 
3  for  the  linear  regressions  with  existing  empirical 
values  as  the  independent  variable,  which  is  expected 
as  they  are  also  functions  of  £.  Selected  regression 
models  that  best  fit  the  data  are  plotted  on  figure  21. 
Regression  [24b]  fits  both  measured  and  §c  and  (J>h  better 
for£<-l  (figs.  21b  and  c),  while  [24c]  appears  to  fit  better 
for  -1  <  £  <  0.  For  cf>m  (fig.  21a),  [24b]  is  plotted  for 
the  west  data  while  [24c]  was  chosen  for  the  east  data. 

Measured  values  of  cj)  for  non-ideal  terrain  can  be 
compared  with  established  values  of  <J>  for  smooth  ter- 
rain using  7  to  illustrate  deviations  in  the  measured 
results.  This  approach  was  previously  established  to 
investigate  the  roughness  sublayer  (Garratt  1978)  but 
also  serves  as  a  useful  comparative  tool  here. 


Ya  = 


fogT  (smooth) 
^qm  (measured) 


[25] 


Figure  22  plots  the  results  for  7m,  7h,  and  7C  vs.  -£  [25]  on 
a  log  scale  using  [16a-b]  for  <j>qT.  The  surface  layer  7h  and 
7m  results  obtained  by  Garratt  (1978)  for  tall  vegetation 
(zD  =  40  cm)  and  by  Fazu  and  Schwerdtfeger  (1989)  for 
bushland  (z0  =  43  cm)  are  also  plotted  on  the  figure  for 
comparison  even  though  the  canopy,  terrain,  and  pro- 
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Table  4.  —  Regression  model  results 


Idl  ulllvlCI  O  iOb 

oorr. 

IVIUUCI 

Wind 

Equation 

dir. 

*, 

b3 

(rho) 

All 

1.8710.12 

7.8610.97 

— 

0.47 

East 

1.9610.18 

6.4111 .23 

— 

0.42 

West 

V  V  CO  I 

1 .7210.17 

11.40+1 .71 

n 

U.JJ 

4>cm  =b,  +  b2K)-b3 

All 

-2.1211.89 

2.5611 .86 

0.0910.05 

0.69 
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-1  50+1  70 

1  95+1  86 

0 1 1 +0  08 

0  88 
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3.2315.01 
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0.66 
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0.4710.05 

0.2410.05 

6.7510.80 

0.71 

East 

0.3910.06 

0.1410.06 

4.3710.72 

0.69 
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V  V  COL 

0  46+0  1 0 

0  39+0  07 

Q 4Q+1  87 

^cm  =b1+exp(b2+b3s) 
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-0.0610.13 

0.6210.08 

2.9310.23 

0.74 

\m  =  b,  +  M 
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1.1310.11 

6.2110.79 

— 

0.47 

East 

1.1110.15 

5.1 010.92 

— 

0.46 
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1  07+0  16 

9  30+1  40 

0  56 

<*>hm  =b1  +  b2K)*a 

All 

0.0410.13 

0.3710.09 

0.37+0.04 

0.74 

East 

-0.4510.27 

0.7610.23 

0.3010.05 

0.79 

West 

0.5110.11 

0.0710.03 

0.6210.08 

0.76 

<l>hm  =  b1  +  exp(b2+b3£) 

All 

0.0610.03 

0.6210.06 

19.012.22 

0.67 

East 

0  55+0  04 

0  78+0  06 

14.211 .45 
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West 

0.8110.05 

1.3810.10 

132.0119.5 

0.72 

<t>hm  =b1  +  exp(b2+b30 

All 

0.1010.13 

0.2610.11 

2.6410.30 

0.62 

<t>mm  =b1  +  b2^ 

All 

0.8910.08 

7.1110.59 

— 

0.82 

East 

1.7910.07 

5.2910.47 

— 

0.71 

West 

0  51+0  06 

4  46+0  61 

0.56 

6     =  b  +  b  (-LYb3 

All 

0  96+0  59 

-0  05+0  54 

0 18+1  33 

0.04 

East 

2.31+3.00 

-1 .3612.97 

-0.0910.23 

0.30 

West 

0.9710.83 

-0.7810.78 

-0.1310.20 

0.43 

<i>mm  =b1  +  exp(b2+b35) 

All 

0.7010.58 

-1.8213.42 

-0.5411.18 

0.08 

East 

1 .0310.03 

-0.6710.18 

1 1 .6014.07 

0.33 
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0.3110.03 
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0.47 

4>mm  =b1  +  exp(b2+b30 

All 

0.7110.05 

-0.7310.14 
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0.33 

file  conditions  at  the  Pawnee  site  are  different.  Figure 
22a  for  ym  illustrates  the  dramatic  effect  of  the  westerly 
winds  on  the  Pawnee  site  (}>m  measurements;  figures  22b 
and  c  do  not  show  this  effect.  With  the  exception  of  west 
7m  values,  the  7  values  are  generally  less  than  1  for  £  < 
0,  indicating  that  in  general  <J>  >  4>qT  for  momentum 
(east  direction),  heat,  and  ozone.  The  7  results  for  £  >  0 
(fig.  23)  are  similar  in  distribution  to  unstable  7C  and  east 
7m  values.  Stable  west  7m  values  (fig.  23a)  are  also 
similar  to  unstable  west  7m  values.  The  stable  7h  results 
are  closer  to  the  ratio  1  but  they  are  widely  distributed. 

If  the  geometric  mean  height,  4.9  m,  selected  as  the 
representative  height  applicable  for  combining  data 
from  3  and  8  meters  were  replaced  with  a  slightly  lower 
value,  say  4  m,  the  7  values  in  figures  22  and  23  would 
increase  and  be  closer  to  the  ratio  1 .  The  choice  of  z  does 


not,  however,  affect  the  eddy-diffusivity  ratio  results 
presented  next. 


Eddy  Diffusivities  of  Ozone,  Sensible  Heat, 
and  Momentum 

Theoretically,  eddy  diffusivities  are  directly  related 
to  atmospheric  stability  and  wind  gradients  and  can  be 
estimated  knowing  u*  without  flux  information.  This  is 
seen  by  combining  [15a]  with  either  [8a],  [8b],  or  [8c]:  Kq 
can  be  calculated  for  any  z,  at  any  stability  £,  knowing  u^ 
and  <pq: 

k  u*  z 

Kq  =  ~ TTrT  [26] 
4  <MsJ 
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Figure  22a — Ratio  of  dimensionless  wind  shear  (momentum)  for 
empirically-established  smooth-terrain  conditions  and  measured 
dimensionless  windshear  [-ym=  <}>m(smooth)/4>m(measured)]  vs. 
stability  parameter  (z/L)  for  unstable  conditions. 
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Figure  22b — Ratio  of  dimensionless  sensible-heat  gradient  for  em- 
pirically-established smooth-terrain  conditions  and  measured 
dimensionless  sensible-heat  gradient  [-yh=  6h(smooth)/ 
4>h(measured)]  vs.  stability  parameter  (z/L)  for  unstable  condi- 
tions. 
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Figure  23a — Ratio  of  dimensionless  wind  shear  (momentum)  for 
empirically-established  smooth-terrain  conditions  and,  measured 
dimensionless  wind  shear  [-ym=  4>m(smooth)Aj)m(measured)]  vs. 
stability  parameter  (z/L)  for  stable  conditions. 
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Figure  23b — Ratio  of  dimensionless  sensible-heat  gradient  for  em- 
pirically-established smooth-terrain  conditions  and  measured 
dimensionless  sensible-heat  gradient  [yh=  <t>m( smooth)/ 
4>  (measured)]  vs.  stability  parameter  (z/L)  for  stable  conditions. 
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Figure  22c — Ratio  of  dimensionless  sensible-heat  gradient  for  em- 
pirically-established smooth-terrain  conditions  and  measured 
dimensionless  ozone  gradient  [-yc=  <t>h(smooth)M>c(measured)]  vs. 
stability  parameter  (z/L)  for  unstable  conditions0. 
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Figure  23c — Ratio  of  dimensionless  sensible-heat  gradient  for  em- 
pirically-established smooth-terrain  conditions  and  measured 
dimensionless  ozone  gradient  [y  =  4>h(smooth)A}>c(measured)]  vs. 
stability  parameter  (z/L)  for  stable  conditions. 


Figure  24  shows  the  resulting  Kq  values  for  (a) 
ozone,  (b)  sensible  heat,  and  (c)  momentum  vs.  £  as 
determined  strictly  from  field  measurements,  [21],  and 
as  determined  using  measured  u*  and  empirical  <J>q  (£), 
[16a-b]  and  [26].  Except  for  measured  K    (west)  and  a 


few  scattered  points,  the  measured  values  of  Kq  tend  to 
be  less  than  the  corresponding  calculated  values.  Again 
the  choice  of  z  plays  a  role  in  calculated  [26].  The 
data  distribution  patterns  vs.  stability  are  very  similar, 
however. 
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Figure  24a — Ozone  eddy-diff usivity  ( Kc)  results;  and  calculated  ozone 
eddy  diffusivity  based  on  measured  friction  velocity  (u.)  vs. 
stability  parameter  (z/L) — May  16-July  11, 1989. 
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Figure  24b— Sensible-heat  eddy-diff  usivity  (Kh)  results;  and  calcu- 
lated sensible-heat  eddy  diffusivity  based  on  measured  friction 
velocity  (u.)  vs.  stability  parameter  (z/L)-May  16-July  11, 1989. 
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Figure  24c— Momentum  eddy-diffusivity  (KJ  results  for  east-  and 
west-wind  categories;  and  calculated  momentum  eddy  diffusiv- 
ity vs.  stability  parameter  (z/L) — May  16-July  11, 1989. 


Calculated  unstable  Kc  andKh  [26]  tend  to  converge 
to  about  20,000  cm2s"1  with  increased  while  mea- 
sured values  tend  to  be  lower  and  scattered  ranging  from 
below  10,000  to  30,000  cm2s  1  for  both,  with  a  few 
outliers  above  30,000  for  measured  Kc.  Stable  eddy- 
diffusivity  values  for  measured  Kc  and  measured  Kh 
follow  the  calculated  values  with  considerably  less  data 


scatter.  Measured  stable  Kc  values  are  generally  less 
than  calculated,  and  measured  stable  values  are 
about  the  same  but  exhibit  more  scatter. 

Calculated  unstable  Km  values  tend  to  converge  to 
about  9,000  cm2s_1  with  increased  while  measured 
values  for  the  east-wind  category  tend  toward  lower 
values:  about  3,000  cm2s"1.  Unstable,  west-wind  cat- 
egory, measured  values  range  from  10,000  cm2s"1  to 
as  high  as  135,000  (there  are  eight  values  beyond 
60,000  cm2s"a  not  shown  on  fig.  24c).  The  swale  west  of 
the  site  is  responsible  for  these  extremely  high  values — 
measured  Au  values  are  considerably  lower  for  the 
west-wind  category  measurements,  resulting  in  higher 
Km  values.  Stable,  measured,  east-wind  category 
values  behaved  similar  to  those  of  Kh  when  com- 
pared to  calculated  values,  while  west-wind  category 
stable  values  are  considerably  higher. 

Table  5  gives  the  linear  regression  results  compar- 
ing measured  Kqm  (y-axis)  and  calculated  KqT  (x-axis) 
values  presented  in  figure  24  (the  second  subscript,  m, 
refers  to  measured  and  T  refers  to  empirical  theory). 
Correlation  coefficients  (table  5)  are  reasonably  good  for 
all  Kh  values  and  stable  values  of  Kc  and  Km.  Surpris- 
ingly, unstable  measured  Kcm  values  for  ozone  did  not 
correlate  well  with  calculated  KcT  values.  Remember 
unstable  cf>c  values  had  better  correlation  coefficients 
compared  to  empirical  theory  than  4>h  values  (see 
table  4). 

The  results  presented  in  figure  24  and  table  5  are  not 
conducive  to  generating  conclusions  concerning  pos- 
sible relationships  between  Kc  and  Km  or  Kh.  Therefore, 


the  ratios  Kc/Km, 


Kh/Km,  and  Kh/Kc,  along  with  the 


empirical  relationships  plotted  in  figure  1,  and  the 
results  determined  from  the  regression  [24a-d]  are  plot- 
ted in  figures  25-27  (note  Kh/Kc  =  (f>c/<J>h,  etc.).  [The 
Kh/Kc  results  are  presented  in  figure  27  for  all  wind 
directions,  the  east-wind  category,  and  the  west  wind 
category,  respectively.] 

The  modeled  eddy-diffusivity  ratio  lines  in  figures 
25-27  are  obtained  by  fitting  the  appropriate  cj>  regres- 
sion, [24a-d],  with  the  related  regression  parameters 
given  in  table  3.  For  instance,  in  figure  25  the  Kc/Km 
ratio  line  (large  dashes)  is  <j>m[24c]/4>c  [24b]  using  the 
parameters  from  table  3  for  the  east  regression  results: 
for  4>m,[24c],  b1  =  1.03,  b2  =  -0.67,  and  b3  =  11.6;  for  cj>c, 
[24b],  ba  =  -1.50,  b2  =  1.95,  and  b3  =  0.11. 

In  figure  25,  the  empirical  ratio  Y^/Km,  (Pruitt  et  al. 
1971)  for  unstable  conditions  doesn't  fit  either  the  east- 
or  the  west-category  data.  The  empirical  ratio  Kh/Km 
(Swinbank  1968)  is  actually  closer  to  the  east  Kc/Km 
regression-based  line  for  unstable  conditions.  The  east 
Kc/Km  regression  line  for  stable  conditions  is  practically 
congruent  with  empirical  Ky/Km.  For  the  west-wind 
category,  the  results  are  dominated  by  large-terrain- 
influenced       values  —  Kc/Km  values  are  all  small. 
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Table  5.  — 


Regression  results:  eddy  diffusivities  of  ozone,  sensible  heat,  and  momentum  vs.  calculated  values. 


KqT  No  Corr. 

V^jnd*    data  coeff. 

Parameter      Stability1-  dir.  slope  SE  const.  SE  points  (rho) 


K    vs  K  T 

cm  cT 


All 
All 
All 
All 


0.37 
0.59 
0.66 
0.60 


±0.06 
0.02 
0.04 
0.03 


3590 
0 

-593 
0 


±5200 
5260 
3410 
3420 


453 
453 
231 
231 


0.26 
0.20 
0.69 
0.69 


Khm  vs  KhT 


All 
All 
All 
All 


0.97 
0.67 
0.74 
0.78 


0.04 
0.01 
0.03 
0.02 


-4920 
0 

402 
0 


3070 
3280 
1990 
2000 


491 
491 
214 
214 


0.76 
0.72 
0.87 
0.87 


K     vs  K 


West 
West 
East 
East 
West 
West 
East 
East 


1.81 
2.54 
0.92 
0.72 
2.33 
2.35 
0.55 
0.63 


0.49 
0.14 
0.10 
0.02 
0.18 
0.10 
0.03 
0.02 


8800 
0 

-1790 
0 

226 
0 

386 
0 


19800 
19900 
3600 
3610 
8860 
8830 
750 
768 


159 
159 
319 
319 
120 
120 
131 
131 


0.28 
0.26 
0.44 
0.44 
0.77 
0.77 
0.82 
0.81 


+  u:  unstable;  s:  stable 

*  East:  030°  -  250°  azimuth  (clockwise) 
West:  250°-030°  azimuth  (clockwise) 
All:     all  wind  directions  combined. 
SE:     Standard  Error 
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Figure  25 — Ratio  of  the  eddy  diffusivities  for  ozone  with  the  east-  and 
west-wind  category  momentum  eddy  diffusivities  (Kc/Km);  theo- 
retical ratio  of  eddy  diffusivities  for  sensible  heat  and  momentum 
(Kh/Km);  theoretical  ratio  of  eddy  diffusivity  for  water  vapor  and 
momentum  (K/KJ;  and  the  ratio  between  regression-model 
east-  and  west-wind  category  dimensionless  wind  shear  and  the 
dimensionless  ozone  gradient  vs.  stability  parameter  (z/L). 


10  - 

9  - 

8  - 

■e- 

7  - 

~~E 

6  - 

5  - 

o 

E 

4  - 

sz 

3  - 

2  - 

1  - 

0  - 

X 

X 

X  ° 
S 

X 

X 

^X 

X  o 

\  o 

o  ^ 

o      °  o 
o       °     "  -  °  o  *>°     •  oo 

<£<*><> 

  i  * — * — IqlrY^H 

  HWI 1  »  . 

-2.6 


-2.2 


-1.8 


-1.4 


-1 

z/L 


-0.6 


-0.2 


0.2 


0.6 


Figure  26 — Ratio  of  the  eddy  diffusivity  for  sensible  heat  with  the 
east-  and  west-wind  category  momentum  eddy  diffusivities 
( Kh/Km);  theoretical  ratio  of  eddy  diffusivity  for  sensible  heat  and 
momentum  (Kh/Km);  and  the  ratio  between  regression-model 
east-  and  west-wind  category  dimensionless  wind  shear  and  the 
dimensionless  sensible-heat  gradient  vs.  stability  parameter 
(z/L). 


(See  Table  4.  for  equation  and  regression  parameters.) 

  Kv/Km(Pruittetal.  1971) 

  Kh/Km(Swinbank  1968) 

  0m  [24b]  /  0c  [24b];  z/L  <  0;  west 

  0m  [24a]  /  0c  [24a];  z/L  >  0;  west 

  0  m [24c]  /  0  c  [24b];  z/L  <  0;  east 

  0m  [24a]  /  0c  [24a];  z/L  >  0;  east 

*  WD  :  West 

O  WD  :  East 

Measurement  error  (precision) 


(See  Table  4  for  equation  and  regression  parameters.) 

  Kh  /  Km  (Businger  et  al.  1971) 

•   Kh/Km  (Swinbank  1968) 

  0m  [24b]  /  0h  [24b];  z/L  <  0;  west 

  0m  [24a]  /  0h  [24a];  z/L  >  0;  west 

  0  m  [24c]  /  0  h[24b];  z/L  <  0;  east 

  0m  [24a]  /  0h[24a];  z/L  >  0;  east 

*  WD  :  West 

O  WD  :  East 

Measurement  error  (precision) 
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(See  Table  4  for  equation  and  regression  parameters.) 

  (Kh  /  Kv  );  see  Figure  1 

  0c  [24a]  /  0h  [24a];  z/L  >  0 

  0c  [24b]  /  0h  [24b];  z/L  <  0 

  0C  [24c]  /  0  h[24c];  z/L  <  0 

  0c[24d]  /  0h[24d];  all  z/L 

*  WD:  West 

O  WD:  East 

Measurement  error  (precision) 
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Figure  27a — Ratio  of  the  eddy  diffusivities  for  sensible  heat  with 
ozone  eddy  diffusivities  (Kh/Kc);  theoretical  ratio  of  eddy  diffusiv- 
ity  for  sensible  heat  and  water  vapor  (Kh/Kv);  and  the  ratio 
between  regression-model  dimensionless  ozone  gradient  with 
the  dimensionless  sensible-heat  gradient  for  all  wind  directions 
vs.  stability  parameter  (z/L). 
10 


o 

b 

o 


9  - 
8  - 
7  - 
6  - 
5  - 
4 

3  H 
2 
1 
0 


(See  table  4.  for  equation  and  regression  parameters.) 

  (Kh/Kv  );  see  Figure  1 

  0C  [24a]  /  0h  [24a];  z/L  >  0;  east 

  <t>c  [24b]  /  0h  [24b];  z/L  <  0;  east 

  0c  [24c]  /  0h  [24c];  z/L  <  0;  east 

O        WD  :  East 


(b) 


-2.6 


-2.2 


-1.4 


z/L 


Figure  27b — Ratio  of  the  eddy  diffusivities  for  sensible  heat  with 
ozone  eddy  diffusivities  (Kh/Kc);  theoretical  ratio  of  eddy  diff  usiv- 
ity  for  sensible  heat  and  water  vapor  (Kh/Kv);  and  the  ratio 
between  regression-model  dimensionless  ozone  gradient  with 
the  dimensionless  sensible-heat  gradient  for  the  east-wind  cat- 
egory vs.  stability  parameter  (z/L). 
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(See  table  4.  for  equation  and  regression  parameters.) 

  (Kh  /Kv  );  see  Figure  1 

  0c  [24a]  /  0h  [24a];  z/L  >  0;  west 

  0c  [24b]  /  0h  [24b];  z/L  <  0;  west 

  0C  [24c]  /  0h  [24c];  z/L  <  0;  west 

*  WD  :  West 


(C) 
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-1.4 
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Figure  27c — Ratio  of  the  eddy  diffusivities  for  sensible  heat  with  the 
ozone  eddy  diffusivities  (Kh/Kc);  theoretical  ratio  of  eddy  diff  usiv- 
ity  for  sensible  heat  and  water  vapor  (K,/Kv);  and  the  ratio 
between  regression-model  dimensionless  ozone  gradient  with 
the  dimensionless  sensible-heat  gradient  for  the  west-wind  cat- 
egory vs.  stability  parameter  (z/L). 

The  Kh/Km  results  in  figure  26  are  similar  to  the 
Kc/Km,  results  except  that  the  east  regression-based 
Kh/Km  line  for  unstable  conditions  increases  dramati- 
cally as  -£  gets  larger.  The  shape  and  location  of  this  line, 
particularly  at  larger  values  of     is  strongly  influenced 


Kh/Km 


by  the  few  measured  4>  values  available  for  larger  The 
ratio  K^/K^  is  expected  to  increase  to  a  constant  value 
as  £  «  -1  based  on  an  argument  that  free  convection 
dominates  under  these  conditions  and  that  turbulence 
is  no  longer  dependent  onu^  (Monin  and  Yaglom  1971). 
Monin  and  Yaglom  (1971)  also  point  out  the  idealized 
possibility  that  K^/K^  -^~as  — >  0  in  free  convec- 
tion. However,  they  quickly  recognize  that  vertical  tur- 
bulence from  thermals  will  continue  to  contribute  to 
K^,  so  that  Kh/Km  should  approach  a  constant  as  £  «  - 
1.  The  results  shown  in  figure  26  are  inconclusive 
relative  to  an  asymptotic  value  for  Kh/Km  as  measured 
£  only  extends  to  -2.6;  however,  the  tendency  for  greater 
values  for  larger  -t,  is  evident  in  the  data. 
The  effect  of  free  convection  has  been  observed  to 
begin  at  £  <  -0.03  (Fazu  and  Scwerdtfeger  1989,  Monin 
and  Yalgom  1971).  In  figure  26  there  appears  to  be  an 
explosion  of  east-category  Kh/Km  values  for  £  <  -0.1, 
which  is  not  evident  for  Kc/Km  (fig.  25).  Values  of  w'd' 
measured  at  the  Pawnee  site  increase  rapidly  at  £  <  -0.1, 
also  marking  the  beginning  of  stronger  heat  fluxes.  For 
stable  conditions  the  Kh/Km  results  are  similar  to  K^/K^ 
with  the  east-category  ratios  slightly  greater  (fig.  26). 

The  Kh/Kc  results  in  three  figures  (figs.  27a-c)  show 
the  effect  of  wind  direction  on  the  regression  results. 
The  ratios  of  the  regression  models  for  the  all-wind 
category  plotted  in  figure  27a  show  a  tendency  for 
Kh/Kc  to  decrease  with  slightly  larger  With  the 
exception  of  the  east-category,  [24b],  regression  ratio  in 
figure  27b,  all  the  data-fitted  regression  ratios  show  the 
same  tendency. This  is  not  immediately  evident  upon 
examining  the  Kh/Kc  data  distribution  by  itself.  If  Kh/Kc 
does  diminish  with  larger  it  is  possible  that  the 
superadiabatic  lapse  rates  associated  with  free  convec- 
tion (strong  solar  insulation  and  low  wind  speed),  could 
be  the  influencing  factor  as  Kh/Kc  is  inversely  propor- 
tional to  Ad.  For  stable  situations,  a  good  average  for 
Kh/Kc  appears  to  be  at  -1.5  for  east-wind  conditions 
and  -1.4  for  west-wind  conditions. 

For  neutral  conditions  as  £  — >  0,  <J)  (£)  should 
approach  the  value  1  (Monin  and  Yalgom  1971).  Hence, 
the  eddy-diffusivity  ratios  should  also  approach  1.  Dif- 
fusivity  ratios  were  tested  for  this  behavior  by  taking  the 
average  ratio  of  all  values  with  I  £  I  <  0.03;  the  numbers 
in  parenthesis  are  the  number  of  cases  averaged  and  the 
±  values  are  standard  deviations  of  the  averages: 


Ratio 

(Kc/Km) 

(Kh  /  Km) 
(Kh/Kc) 


Wind  Direction  Category 
All  East  West 

0.64  ±  cp^  53  1.16  ±  0.53  0.33  ±0.28 
(118)             (34)  (79) 

0.66  ±0.48  1.15  ±  0.56  0.46  ±  0.25 
(103)             (29)  (74) 

1.34  ±  0.76  1.07  ±  0.56  1.45  ±  0.80 
(91)              (27)  (64) 
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The  east-category  neutral-stability  results  are  very  close 
to  the  value  1  for  all  three  ratios,  while  the  west-category 
results  involving  Km  are  considerably  less  than  1, 
showing  the  effect  of  small  measured  Au  values. 

The  west-category  neutral-stability  result  for  Kh/Kc 
is  considerably  larger  than  the  east-category,  which 
appears  to  contradict  the  claim  that  the  measured  values 
involved  are  not  wind-direction  dependent.  Examina- 
tion of  the  individual  records  used  for  these  calcula- 
tions revealed  that  there  are  a  larger  number  of  small 
positive  Ad  values  as  a  result  of  editing  for  the  west- 
category  compared  to  the  east-category,  thereby  artifi- 
cially raising  the  average.  Figure  27a  also  shows  that 
Kh/Kc  values  for  the  west- wind  category  close  to  £  =  0  are 
somewhat  larger  but  as  a  whole  are  distributed  similar 
to  the  east-wind  category. 

The  plotted  Kh/Kc  values  in  figure  27a  appear  to 
take  a  step  jump  at  £  =  0,  with  higher  values  for  £  >  0  and 
lower  values  for  £  <  0.  The  Kh/Kc  values  also  tend  to 
converge  from  a  broad  spread  of  data  at  £  =  0  to  a  value 
-1.5  as  £  — »0.4,  the  largest  £  values  in  this  analysis.  This 
behavior  is  not  as  evident  for  Kc/Km  (fig.  25)  or  for  Kh/ 

(fig.  26).  The  effect  is  controlled  by  Ac/A0;  Ad  is 
very  small  relative  to  Ac  near  £  =  0,  then  A  6  gradually 
increases  to  provide  a  relatively  fixed  ratio  for  Ac/A0 
between  2  and  3  ppb/°C  for  £  >  0.2. 


CONCLUSIONS 

The  reliability  and  validity  of  the  Pawnee  site  flux 
and  gradient  data  are  established  in  the  previous  sec- 
tions and  in  Zeller  et  al.  (1989).  In  some  cases  the 
Pawnee  site  results  verify  existing  knowledge,  and  in 
other  cases  the  results  are  site-specific  and  deviate  from 
the  results  of  other  investigators.  The  conclusions  of 
this  study  relative  to  the  Pawnee  site  measurements  and 
to  the  vertical  turbulent  transport  of  ozone  are  given 
below.  Table  6  provides  a  summary  of  the  eddy-diffu- 
sivity  ratio  results  for  Kh/Kc,  Kh/Km  and  Kc/Km. 


Eddy  Diffusivity  Comparisions 

The  data  and  regression  model  eddy-diffusivity 
ratio  results  presented  in  figures  25-27,  and  the  dimen- 
sionless  gradient  and  regression  results  presented  in 
figure  21,  provide  the  basis  for  conclusions  (table  6).  The 
values  of  c}>q  for  ozone,  sensible  heat,  and  momentum 
(even  though  the  momentum  results  are  wind-direction 
sensitive),  are  well  behaved  and  add  credence  to  the 
broadly  distributed  eddy-diffusivity  ratios  presented  in 
figures  25-27.  The  (J>  values  obtained  using  data  from 
two  heights  are  sensitive  to  choice  of  z  =  4.9  m,  which 
is  evident  in  the  plotted  values  of  ^q  (figs.  22-23).  It  is 
also  evident  that  Monin-Obukhov  (M-O)  similarity  does 
not  completely  explain  the  behavior  of  <j>  for  west 
winds. 

For  Kh/Kc  (fig.  27a)  the  conclusion  Kc=  Kh  for 
unstable  cases  is  reasonable  within  the  experimental 
error  of  the  study.  The  regression  models  for  Kh/K 
(=<j>c/<J)jl)  call  attention  to  the  possibility  that  Kh/Kc 
actually  decreases  with  This  is  physically  plausible 
if  w'c' ,  w'  6'  and  A  c  are  limited  while  A  6  increases  with 
free  convection.  Since  Ad  is  also  eventually  limited,  a 
lower  asymptotic  limit  with  £«-l  is  most  probable 
however. 

For  the  stable  case,  results  indicate  Kc*  Kh.  The 
appropriate  conclusion  is  that  on  the  average  Kh—  1 . 5  Kc 
for  east  winds  and  Kh—  1.4  Kc  for  west  winds.  It  is  most 
likely  for  stable  situations  that  the  chemical  reactivity  of 
ozone  with  surface  constituents  is  playing  a  role  and 
that  this  result  is  specific  to  ozone.  During  stable  periods 
ozone  is  not  replenished  from  above  as  rapidly  as  during 
unstable  periods.  This  allows  the  ozone  in  the  lowest 
layer  to  be  depleted  chemically  while  this  would  not 
happen  to  a  nonreactive  mass,  hence  larger  Ac  values 
for  ozone  and  consequently  lower  Kc  values;  therefore 
larger  Kh/Kc  ratios.  The  waning  in  Kh/Kc  data  scatter  as 
£  increases  from  0  supports  the  above  argument  as  it 
indicates  that  a  transition  is  taking  place  for  measured 
A0/Ac  as  stability  intensifies. 


Table  6.  —  Eddy  diffusivities  ratio  results. 


Ratio 

Unstable 

Stable 

East* 

West 

East 

West 

(smooth 

(irregular 

(smooth 

(irregular 

terrain) 

terrain) 

terrain) 

terrain) 

n  c 

K.  -  1 .5K 

n  c 

Kh=1.4Kc 

Kh/K 

h  m 

+ 

Kh  -  1 .35K 

h  m 

Kr/Km 

c  m 

+ 

Kc=0.35Km 

K  -  0.9K 

c  m 

K-0.35Km 

*  East:  wind  direction  030°  -  250°  Azimuth. 
West:  Wind  direction  250°-030°  Azimuth. 
+  conclusions  not  possible. 
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The  regression  results  for  Kk/Km  =  a  (east-wind 
direction),  as  shown  in  figure  26  for  unstable  cases, 
provides  an  interesting  argument  for  the  theoretical 
free-convection  ratio  Kh/Km  — >  ©°.  There  are  too  few  data 
points  to  take  this  conclusion  seriously:  however,  mea- 
sured Kh/Km  are  definitely  larger  than  results  presented 
for  smooth  terrain.  A  plausible  explanation  is  that  Au 
values  are  affected  (increased)  by  the  -1%  terrain  slope 
experienced  by  east  winds.  For  stable  Kt/Km  (east-wind 
direction)  the  measured  results  plotted  in  figure  26 
indicate  an  average  of  1.3:  Kh/Km  ->  1.3  as  £  ->  0.4. 
Although  the  data  scatter  is  significant,  it  is  a  reasonable 
result  for  stable  cases  where  the  constant  a  ~  1.35  was 
expected.  The  results  for  K^/K^  ~  a  (west-wind  direc- 
tion) are  significant  in  that  they  demonstrate  that  the 
effect  of  irregular  terrain  is  to  lower  the  ratio  of  K^/K^. 
This  is  plausible  as  would  increase  with  increased 
turbulence  mechanically  generated  by  irregular  terrain 
features.  The  average  result  K^/K^ 5=3  0.5,  a  constant  for 
both  stable  and  unstable  conditions,  should  be  restricted 
to  the  lower  surface  layer  as  the  irregular-terrain  effect 
would  lessen  with  increased  z. 

The  conclusions  for  Kc/Km  (fig.  25)  are  similar  to  the 
conclusions  for  Kh/Km.  The  data  distribution  shown  in 
figure  25  is  very  similar  to  figure  26,  although  the  data 
values  are  different  enough  to  generate  a  plausible 
regression-based  Kc/Km  result  for  the  unstable  east- 
wind  category  (i.e.,  Kc/Km  approaches  an  asymptotic 
limit  as  £  — >  -«>).  The  regression-based  Kc/Km  in  figure  25 
for  I  <  0  isn't  too  different  from  the  Kh/Km  ratio  provided 
by  Swinbank  (1968).  For  the  stable  east-wind  category, 
the  average  KJK^  —  0.9,  which  is  very  close  to  the  result 
provided  by  Pruitt  et  al.  (1971)  for  K/K^.  For  both 
stable  and  unstable  west-wind  categories,  the  average 
ratio  Kc/Km  —  0.35  is  perhaps  a  useful  result  for  the 
lower  surface  layer  in  irregular  terrain. 

The  comparison  Kc  =  K^,  with  (different  from  Kh) 
based  on  the  empirical  expression  given  in  [11]  seems 
reasonable  for  stable  conditions  but  not  for  unstable 
conditions  (with  the  exception  of  the  east-wind  regres- 
sion model  [24b])  based  on  the  results  shown  in  figures 
25  and  27a.  Because  of  the  chemical  nature  of  ozone  this 
conclusion  is  specific  to  ozone. 

Table  6  provides  a  summary  of  the  conclusions 
discussed  above.  All  results  for  stable  conditions  are 
consistent  within  respective  wind-direction  categories 
(i.e.,  Kh/Kc  =  (Kk/KJ/tKc/KJ).  There  are  no  conclu- 
sions for  the  unstable  east  category.  The  unstable  west 
category  conclusions  are  not  consistent  within  them- 
selves (i.e.,  Kh  *  Kc  if  Kh  -  0.5Km  and  Kc  -  0.35KJ 
because  the  results  Kh  —  Kc  for  unstable  west-wind 
conditions  are  based  on  data  that  for  the  most  part  are 
associated  with  near  neutral  stability,  £  >  -  0.2  (fig.  27c); 
they  are,  however,  the  best  conclusions  that  can  be 
reached  given  the  distribution  in  eddy-diffusivity  ratio 
results. 


Constant  Flux  Layer 

The  8-  and  3-meter  measurement  heights  used  at  the 
Pawnee  site  are  too  closely  spaced  to  expect  detection  of 
any  significant  flux  differences;  hence,  the  terms  "con- 
stant flux"  and  "constant  flux  divergence"  are  synony- 
mous for  this  analysis.  The  measurement  results  for 
corrected  fluxes  (fig.  20)  compared  to  those  for 
uncorrected  fluxes  (fig.  16),  show  that  the  constant-flux 
layer  assumption  for  sensible  heat  and  momentum  is 
reasonable,  as  the  results  are  essentially  the  same  for 
both. 

For  ozone,  the  constant  flux  or  flux-divergence  layer 
assumption  is  not  as  easy  to  evaluate.  Figure  17,  the 
results  for  uncorrected  ozone  flux  (corrected  for  heat 
flux),  shows  that  there  is  definitely  less  flux  measured 
at  the  3-meter  height  compared  to  the  8-meter  height. 
These  data  show  strong  flux  divergence  in  the  down- 
ward direction  (i.e.  less  flux  at  3  meters  compared  to  8) 
as  would  be  expected.  This  is  based  on  assuming  ozone 
is  supplied  from  above  8  meters  and  is  not  generated 
between  8  and  3  meters .  It  is  unknown  in  this  study  if  the 
measured  rate  of  flux  divergence  for  w'c'  would  con- 
tinue above  8  meters;  however,  it  is  unlikely  based  on 
aircraft  measurements  at  higher  elevations  that  report 
flux-related  values  for  ozone  in  the  same  range  as  mea- 
sured at  the  Pawnee  site  (Lenschow  1982).  The  results 
for  corrected  ozone  flux  shown  in  figure  18,  obtained  by 
applying  all  known  and/or  theorized  flux  corrections, 
show  that  a  constant  flux  layer  (as  opposed  to  a  constant 
flux-divergence  layer)  for  ozone  is  statistically  accept- 
able in  that  w'c' (8  m)  =  w'c' (3  m)  within  a  95%  confi- 
dence interval. 

The  results  in  figure  18  are  quite  scattered  and  are 
produced  by  flux  corrections  that  average  +45%  for  8 
meters  and  +85%  for  3  meters.  The  corrections  for  stable 
cases  are  particularly  large  and  uncertain,  and  the  cor- 
rections for  unstable  cases  are  large  compared  to  the 
heat-flux  correction  for  w'c'  (maximum:  -25%).  After 
evaluating  the  corrected  and  uncorrected  ozone  flux, 
there  are  still  two  questions:  how  valid  are  the  flux 
corrections;  and,  can  reactions  with  NOx  species,  aero- 
sols, and  NMHC's  above  3  meters  really  be  discounted? 
For  these  reasons  the  constant  flux  or  the  constant  flux- 
divergence  question  for  ozone  flux  is  still  open  to 
further  research.  However,  assuming  the  flux  correc- 
tions tested  are  at  least  partially  correct  and  that  net 
ozone  production  between  3  and  8  meters  is  minimal, 
the  assumption  of  a  constant  flux  for  ozone  is  reason- 
able. 


Other  Findings 

The  vertical  transport  and  gradients  of  ozone  and 
sensible  heat  behave  similarly.  Dimensionless  gradient 
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values,  4>  ,  for  both  parameters  follow  M-O  similarity. 
The  gradient  of  ozone  has  a  diurnal  pattern  that  appears 
to  be  controlled  by  a  combination  of  sources  and  sinks 
as  well  as  boundary-layer  mixing  height. 

Terrain  features  have  a  dramatic  effect  on  Km  and 
4>mby  affecting  Au,  while  terrain  effects  on  z0  (calcu- 
lated with  [17a]  and  t,  as  well  as  Kc,  K^,  <$>c,  and  fy^  are 
not  discernible. 

Flux-measurement  corrections  for  inadequate 
sensor  response  are  measurement-height  sensitive 
and  play  a  factor  in  measured  results  at  lower  measure- 
ment elevations  for  ozone.  Flux  corrections  for  sensible 
heat  flux  and  momentum  flux  are  not  as  sensitive  to 
elevation. 

Recommendations  for  Further  Research 

The  following  list  of  questions,  based  on  the  results 
of  this  investigation,  provide  recommendations  for  fu- 
ture research. 

•  Does  the  ratio  Kc/Kh  for  ozone  and  heat  dimin- 
ish with  increasing  -£? 

•  Is  the  measured  flux  divergence  for  ozone  within 
the  surface  layer  as  strong  as  indicated  by  the 
Pawnee  site  measurements  and  how  high  does  it 
extend  before  becoming  constant?  (The  answer 
will  require  simultaneous  flux  measurements 
from  more  than  two  heights.) 

•  What  controls  the  transition  of  A  c  /A  6  for  0  >  £ 
>  0.2  from  widely  fluctuating  values  to  a  rela- 
tively stable  ratio? 

•  How  important  are  the  response-related  flux 
corrections  for  fluxes  measured  under  stable 
conditions,  and  what  can  be  done  to  improve 
the  corrections  in  general? 

•  Does  Kh/Kc  for  a  nonreactive  depositing  mass 
tracer,  where  c  is  not  ozone,  behave  differently 
from  Kh/Kc  where  c  is  ozone? 
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APPENDIX  A. 
DATA  ACQUISITION  AND  HANDLING 

Data  Acquisition  and  Real-Time  Computations 

Data  collection  and  processing  were  accomplished 
using  a  PC  (Compaq  386)  interfaced  with  an  analog-to- 
digital  conversion  board  (DT2805).  Eight  channels  were 
available  on  each  of  two  boards  for  data  acquisition.  A 
total  of  16  channels  (table  A-l)  were  sampled  at  ca.  14 
Hz.  Output  measures  include  means,  variances,  covari- 
ances,  and  other  micrometeorological  quantities.  The 
software  used  for  real-time  data  processing  was  adapted 
from  the  BASIC  flux  program  for  eddy  correlation  in 
nonsimple  terrain  developed  by  McMillen  (1986).  The 
compiled  flux  program  is  capable  of  running  the  entire 
eddy-correlation  measurement  system  unattended  for 
several  days.  Tasks  accomplished  in  real-time  by  this 
program  include  real-time  data  acquisition  and  process- 
ing, raw  data  storage  if  required,  statistical  summaries  at 
the  end  of  the  assigned  averaging  period  (usually  one- 
half  hour),  and  instrumentation  zero  each  half  hour  to 
ensure  accurate  chemical  measurements.  Fluctuations 
for  each  data  point  (c'(t),  0'(t),  w'(t),  etc.)  were  calcu- 
lated in  real-time  by  subtracting  the  instantaneous  sig- 
nal from  a  200-second  running  mean  generated  by  using 
a  recursive  digital  filter: 

(         .  Y 


c(t)  =  [c(t)/(200*n)] 
c'(t)  =  c(t)-c(t) 


:(t-l) 


200*n  ) 


[A-l] 


where  n  =  sampling  frequency. 

The  value  200  seconds  was  chosen  to  maximize  the 
covariance  based  on  test  calculations  and  experience 
with  other  surface-layer  eddy  correlation  experiments 
(McMillen  1986).  The  covariance  calculations  involv- 
ing wind  and  ozone  concentration  were  handled  in  real- 
time by  saving  the  wind  data,  u(t),  v(t),  and  w(t)  for 
lagtime  *  n  data  collection  loops.  For  example,  if 
lagtime  =  2  seconds  and  the  data  were  collected  at 
14  Hz,  14  *  2  =  28,  then: 


w'c'(t)  =  w'(t-28)c'(t) 
u'c'(t)  =  u'(t-28)c'(t) 
v'c'(t)  =  v'(t-  28)c'(t) 


[A-2] 


where  t  =  1  is  equivalent  to  1/14  second. 

At  the  end  of  each  half-hour  sampling  period,  the 
mean  horizontal  wind  direction  (resulting  v  =  0), 
77;  the  vertical  mean  wind  direction  (mean  streamline, 
resulting  w=  0),  0;  and  the  streamline  "twist"  angle 
(resulting  wV  —  0),  /3,  were  calculated  and  the  sampled 
covariance  matrix,  [M],  was  then  rotated  three  times  (17, 
d  and  (3)  so  that  the  output  covariances  were  aligned 
with  the  mean  wind  streamline  (Wesely  1970). 


[M]  = 


,/2 


w'*  w'u'  w'v'w'T'w'c' 


u'w'  u'2    u'v'  uT'  u'c' 


v'w'  v'u'  v'2  v'T'  v'c' 


[A-3] 


While  the  above  calculations  were  made  and  all 
summary  data  (table  A-l)  were  stored,  the  chemical 
instruments  were  zeroed  and  all  parameters  were 
reinitiated  for  the  next  sample  period.  The  chemical 
instruments  were  then  stabilized  with  ambient  air  prior 
to  commencement  of  the  next  data-collection  period. 
The  above  tasks  required  ca  1.8  minutes  per  half-hour 
data-acquisition  cycle. 


Data  Editing 

Data  handling  and  editing  for  eddy-correlation  mea- 
surements must  be  accomplished  with  great  care.  The 
flux  measurement  is  a  statistical  covariance  that  be- 
comes meaningless  if  either  timing,  calibrations,  or 
calculations  are  off.  In  addition  to  the  quantities  listed 
in  table  A-l,  calibration  values,  calibration  offsets,  and 
the  variance  of  each  covariance  are  also  calculated  and 
stored.  There  were  171  data  averages  archived  every 
half  hour.  Data  handling  must  be  carefully  monitored. 
Prior  to  testing  for  stationarity  and  homogeneity,  to 
eliminate  data  that  may  have  contained  bogus  or  corrupt 
information,  records  that  met  the  criteria  listed  in  table 
A-2  were  eliminated.  Korrell  et  al.  (1982)  used  similar 
wind-related  criteria  for  their  surface-layer  analysis  of 
Boulder  tower  data,  Fazu  and  Schwerdtfeger  (1989)  also 
applied  similar  data-editing  criteria  for  wind  and  tem- 
perature flux  and  gradient  data  collected  over  bushland 
in  Australia. 

Ozone  Gradient  Determination 

A  correct  Ac  measurement  for  ozone  is  important 
for  the  eddy-diffusivity  calculation  [20].  This  analysis 
used  A  c  data  determined  from  two  TECO  ozone  analyz- 
ers sampling  from  two  different  heights,  rather  than  the 
traditional  approach  of  switching  the  inlet  of  a  single 
analyzer  between  two  levels,  which  leads  to  errors 
(Woodruff  1986).  The  Ac  measurement  technique  used 
alternated  the  level  each  TECO  monitored  every  half 
hour  through  the  use  of  two  three-way  computer-con- 
trolled solenoid  switches.  Since  each  CAAM  was 
plumbed  to  a  single  level  and  not  allowed  to  switch 
between  levels  during  any  1-week  sampling  period 
while  the  TECOs  did  switch,  four  sets  of  instrument 
regression  calibrations  similar  to  figure  5  could  be 
extracted:  CAAMl:TECOl;  CAAMl:TEC02;  CAAM2 
:TECOl;  and  CAAM2:TEC02.  Hence,  two  separate  cali- 
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Table  A-1 .  —  Data  acquisition  and  real-time  measured  quantities  at  the  Pawnee  site. 


DT  2805 

channel  Quantity 
no. 

Height  (m) 

Output 
(every  1/2  hour) 

1 

w:  vertical  wind  speed  component 

8 

—  ~/i 
w,w 

2 

u:  horizontal  wind  speed  component 

8 

—  ,2 

u,u 

3 

v:  horizontal  wind  speed  component 

8 

-  ,2 
V,V 

4 

T:  temperature  (6  -  1 .05  T) 

8 

T.f2 

5 

c:  (TEC049):  ozone 

8 

c 

6 

c:  (CAAM):  ozone 

8 

c,c'2 

7 

water  vapor 

8 

q,q'2 

8 

open 

9 

w:  vertical  wind  speed  component 

3 

w,w'2 

10 

u:  horizontal  wind  speed  component 

3 

—  ,2 

u,u 

11 

v:  horizontal  wind  speed  component 

3 

-  ,2 
V,V 

12 

T:  temperature 

3 

13 

c:  (TEC049):  ozone 

3 

c 

14 

c:  (CAAM):  ozone 

3 

c,c'2 

15 

open 

16 

open 

* 

Momentum  flux 

3&8 

w'u'.w'v' 

* 

Heat  flux 

3&8 

w'0' 

* 
* 

Ozone  flux 

Horizontal  wind  speed 

3&8 
3&8 

w'c' 
(U2+v2) 

* 

Friction  velocity 

3&8 

l-w  u  I 

* 

Drag  coefficient 

3&8 

(-w'u'j/  u2 

Horizontal  wind  rotation  angle 

3&8 

* 

Vertical  wind  rotation  angle 

3&8 

e 

* 

Streamline  twist  angle 

3&8 

* 

Sample  frequency 

3&8 

* 

No.  of  samples/sample  period 

3&8 

N 

* 

Covariance  rotation  matrix 

3&8 

[M] 

21X 

Temperature  gradient 

3&8 

Data  calculated  from  data  recorded  for  each  height. 
21 X      Data  acquired  with  Campbell  model  21 X  data  logger. 


bration  or  alignment  could  be  developed  for  TECOl 
compared  to  TEC02: 


TECOl 


CTEC02^AAMl^     anC*      C  TECOl  - 


CTEC02(CAAM2). 

The  results  for  the  week  June  27-July  4,  1989,  are  given 
as  an  example: 


c TECOl  =  °-927  CTEC02(CAAM1)  -1.915  [A-4a] 

CTEC01  =  0.920  CTEC02(CAAM1)  -1.927  [A-4b] 

These  equations  are  practically  congruent  and  thus  both 
verify  the  validity  of,  and  give  the  inherent  error  in  this 
approach  for  aligning  the  two  TECOs;  that  is,  [A-4a]  and 
[A-4b]  are  arrived  at  independently.  The  two  equations 
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Table  A-2.  —  Data  editing  criteria. 


Quantity 


Criteria 


Reynolds  stress 

101:  vertical  angle 

1/81:  twist  angle 

TECO  variance 

CAAM  variance 

CAAM  covariance 

CAAM  covariance  variance 

Temperature  variance 
Temperature  covariance 
Temperature  covariance  variance 

w' 


'6'/(dT/d\ 
vFc'/(dc/dt 


\Ad\ 

Wind  speed* 

Wind  Speed 

Wind  speed  difference 

Wind  direction 

lAul 

Instrument  malfunction 
Instrument  malfunction 
Instrument  malfunction 
Instrument  malfunction 


>  0.0  m2  s2 
>+  10° 

>  +  20° 

>  10  ppb2 

>  0.2  volts2 

8  m  direction  *3m  direction 

>  0.011  (volts  m)2s  2 

>  10.00 

>  1 .2  °C2 

8  m  direction  *  3  m  direction 

>  0.4  (°C  m)2  s"2 

<  z 

<  z 

<0.05 

Raining  or  not 

>  1 5  m  s~1 
u(3)  u  u(8) 
230°  >  WD>  190° 
<0.05 
Electronic 
Eoysn-y  depleted 
Calibration 
Broken  Sensor 


Typical  no.  of 
values/week 

15 

11 

6 

15 

5 

5 

3 
16 

4 
25 

1 

12 

Occasional 
3 

Variable 
Occasional 
20 
10 

3 

Occasional 
Occasional 

2 

Occasional 


*  Sonic  anemometers  work  incorrectly  when  wet.  Wind  speeds  >  15  m  s~1  infer  gusts  >  20m  s_1  , 
the  upper  limit  for  ATI  sonic  anemometers. 


are  next  averaged  and  then  one  of  the  TECOs  is  aligned 
with  the  other;  e.g.,  the  result  from  [A-4a]  and  [A-4b]: 

CTEC02 (aligned)  =  0.923  CTEC02  (original) -1.921  [A-5] 

In  [A-5]  cTEC02  (aligned)  is  the  adjustment  of  cTEC02 
(original)  required  to  bring  it  into  alignment  with  c  TECOl. 
Once  this  procedure  is  accomplished,  Ac  can  be  calcu- 
lated by  subtraction.  An  example  of  the  result  of  this 
procedure  for  May  19th,  1989  (Julian  day  139)  is  shown 
in  figure  A-l.  The  saw-tooth  effect  of  the  Ac  an  indica- 
tion that  the  TECOs  were  still  not  completely  aligned. 
Since  the  TECOs  are  switched  every  other  half  hour,  the 
use  of  the  information  within  the  differences  is  in  order. 
Equation  [A-6]  was  applied  to  the  aligned  Ac  data;  it  is 
equivalent  to  smoothing  over  three  points  with  weights 
1/4,  1/2,  and  1/4: 


The  result  is  the  dashed  line  in  figure  A-l. 

A  separate  4-day  test  was  run  on  the  Ac  ozone  data 
collection  and  calculation  scheme  described  above. 
Commencing  midday  July  22, 1989  (Julian  day  202),  the 
intakes  of  both  TEC049s  were  connected  to  the  same 
intake  manifold;  hence,  the  same  level  (8  meters)  was 
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Figure  A-1 — Detail  of  ozone  gradient  data  and  smoothing  (switch 
compensated)  calculation  results — May  19, 1989. 
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sampled  by  both  TECOs.  In  theory  Ac  should  be  zero 
for  the  period  of  time  both  TECOs  were  sampling  the 
same  air.  Ac  data  prior  to  and  after  midday  July  22, 
1989,  demonstrates  the  accuracy  of  the  overall  Ac 
measurement  scheme:  for  the  4-day  test  period  shown  in 
figure  A-2,  Ac  =  +  0.07  ppb  and  crA-  =  0.13  ppb 
(variance  =  0.017  ppb2) ,  which  represents  the  bias  in  the 
system  and  the  precision  of  the  Ac  measurement 
for  Az  =  5m. 


Flux  Measurement  Corrections 

Although  eddy  correlation  is  the  flux-measurement 
technique  of  choice,  recent  attention  has  been  given  to 
the  requirement  for  careful  data  processing  of  eddy- 
correlation  data  (Businger  1986;  Wesely  et  al.  1989). 
Businger  (1986)  lists  ten  specific  items  related  to  flux 
measurements  that  researchers  in  the  past  have  not 
diligently  addressed.  In  a  plea  for  flux  research  unifor- 
mity, he  proposes  that  the  flux  correction  items  listed  in 
table  A-3  at  least  be  noted  in  future  research  reporting, 
even  if  they  are  not  used,  so  that  results  can  be  evaluated 
on  a  comparative  basis.  As  indicated,  nine  applicable 
areas  of  concern  for  gaseous  eddy-correlation  measure- 
ments were  considered  for  the  Pawnee  site  data  ana- 


lyzed herein.  Except  for  items  1,  2,  6,  7,  8,  and  10a,  the 
concerns  in  table  A-3  are  not  routinely  corrected  for 
because  they  are  theoretical  in  nature  and  add  consider- 
able scatter  to  the  data  (Kaimal  1989;  Wyngaard  1989a). 
Flux  corrections  for  items  3  and  4  involve  instrument 
frequency  response;  these  corrections  are  complicated 
and  are  presented  in  Zeller  et  al.  (1989). 


Effects  of  Heat  and  Water  Vapor 

Webb  etal.  (1980)  provide  corrections  ([24],  Webb  et 
al.  1980)  for  the  effects  of  heat  flux  and  water  vapor  flux 
on  mass-flux  measurements  made  with  instruments 
that  measure  density  fluctuations  as  opposed  to  mixing- 
ratio  fluctuations.  Essentially,  gradients  of  heat  and 
water  vapor  can  cause  differences  in  air  density  moving 
upward  versus  downward.  Because  the  correction  for 
heat  is  about  5  times  that  of  water  (Webb  et  al.  1980)  and 
because  relative  humidity  is  generally  low  at  the  Paw- 
nee site  [A-7],  which  accounts  for  the  heat-flux  correc- 
tion only,  is  used  to  correct  ozone-flux  measurements 
made  with  the  CAAM  instruments  (Zeller  et  al.  1989): 


Fr  =  Fr  (measured )  +  c  ^  ^ 

e 


[A-7] 


Table  A-3.  —  List  of  dry  deposition  measurement  concerns  applied  to  the  Pawnee  site  data  set 

(Businger  1986). 


Flux  Corrections 

Addressed 
(yes/no) 

Affected 
measurement(s) 

1 .     Effects  of  heat  and  water  vapor  § 

yes:  ( w'0'  only) 

w'c' 

2.     Insufficient  averaging  time 

yes+:  28.2  min 

means,  fluxes 

3.     Inadequate  sensor  response: 

a.  dynamic  frequency  response 

yes* 

fluxes,  variance 

b.  sensor  mismatching 

yes* 

c'w' 

c.  line  averaging 

yes* 

u'w',  w'0' 

d.  data  logging 

yes* 

all  instruments 

e.  aliasing 

yes* 

fluxes,  variance 

f.  volume  averaging 

yes* 

w'c' 

4.     Sensor  separation 

yes* 

w'c',  w'0' 

5.     Random  noise 

no 

6.     Homogeneity  and  stationarity 

yes+ 

fluxes 

7.     Sampling  height 

yes+ 

fluxes 

8.     Inadequate  fetch 

yes+ 

all  measures 

9.     Deliquescences  of  particles 

NA 

NA 

10.     Flow  distortion 

a.  mounting  structures 

yes+ 

vector  quantities 

b.  instrument  shadow 

yes* 

vector  quantities 

§  The  heat  flux  portion  of  this  correction  was  applied  to  all  ozone  fluxes:  both  "corrected"  and 

"uncorrected"  fluxes  referred  to  in  the  text  have  been  corrected  for  heat  flux. 
*    Evaluated  but  not  corrected  for  in  the  eddy  diffusivity  analyses. 
+  Part  of  data  editing  criteria  (see  table  A-2). 
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where  6  is  the  average  potential  temperature  in  °K. 

The  Webb-Pearmen-Luening  correction  (Webb  et  al. 
1980),  referred  to  herein  as  the  "heat-flux  correction, "  is 
the  only  one  routinely  made  so  that  Fc  (measured) 
in  [A-7]  is  the  only  flux  measurement  that  is  actually 
adjusted  on  a  routine  basis.  These  corrections  are  gener- 
ally positive  and  between  0%  and  25%  (0  to  0.04  ppb  m 
s"1)  of  the  measured  ozone  flux.  The  resulting  corrected 
values  for  w'c'  are  therefore  less  negative  or  smaller  in 
absolute  value. 

Insufficient  Averaging  Time 

The  appropriate  length  of  time  for  a  good  average  in 
field  turbulence  measurements  is  a  balance  between 
capturing  a  statistically  significant  number  of  important 
eddies  in  the  dominant  scale  while  atmospheric  condi- 
tions remain  stationary  and  homogeneous.  Businger 
(1986)  has  reviewed  this  problem  and  offers  the  averag- 
ing period  scale  T-  as  an  empirically  derived  value  for 
estimating  the  averaging  time  to  determine  c  with  an 
accuracy  of  a: 


T-  = 

1  c 


20z  c/2 

2—  -2 

a  u  c 


[A-8] 


For  T^:,  the  corresponding  averaging  time  to 
determine  w'c'  with  an  accuracy  of  a  (a  =  0% ,  unattain- 
able in  practice,  means  no  error): 

(   A 

20z 


w  c 


a2u 


(w'c'f 


2  2 
U  C 


-1 


[A-9] 
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Figure  A-2 — Demonstration  of  the  accuracy  and  precision  of  the 
ozone  gradient  measurement — Ozone  gradient  test:  July  21-25, 
1989. 


Random  Noise 

Flux  data  gathered  at  the  Pawnee  site  are  not  ad- 
justed for  random  noise.  With  the  exception  of  not 
employing  a  Butterworth  filter,  a  significant  effort  was 
made  to  eliminate  noise  from  all  sensors  and  data  lines. 
Ground  loops  were  avoided  by  grounding  signal  lines 
on  only  one  end,  the  shelter  and  all  systems  were 
grounded,  and  the  data-collection  system,  PC/DT2805, 
was  operated  through  an  uninterruptable  power  supply 
which  also  provided  a  constant  power  source.  All  16 
input  channel  data  were  plotted  on  the  PC  screen  on  a 
continuous  basis.  In  this  way  spurious  signals  and 
excessive  noise  in  the  data  were  visually  detected  and 
eliminated  early  in  the  data-collection  process.  Also, 
the  system  was  monitored  weekly  for  any  changes  in 
sensor  performance  (Stocker  and  Lukens  1988). 


Tincreases  linearly  with  z;  therefore  for  z  =  8  m  and 
T-  =  28.2  minutes,  typical  average  measured  val- 
ues of  c'2  /  c2  —  0.0017,  and  u  -  4.0  m  s"1  for  the  week 
June  13-20,  1989,  as  an  example,  gives  a  ~  0.6%,  the 
accuracy  with  which  c  can  be  measured  for  the  Pawnee 
site.  Using  [A-9]  with  typical  values  u*  =  0.27  m  s"1, 

c*=0.35  ppb,  and  (w'c')2  =  0.00016  m2  ppb2  s"2  for  the 

same  week  in  June  gives  a  ~  13.7%,  the  accuracy  with 
which  w'c'  can  be  expected  to  be  measured  within  the 
same  28.2-minute  averaging  period.  The  T^p-,  value  for 
a,  13.7%,  is  very  reasonable,  as  Businger  (1986)  consid- 
ers values  of  a  =  100%  for  T—r~,  workable  based  on  the 

_     w  c 

fact  that  the  accuracy  of  c  measurements  quickly  dete- 
riorate when  attempting  to  extract  a  gradient,  Ac/Az. 
This  expectation  is  true  for  the  Pawnee  site,  based  on  the 
data  presented  in  figure  A-2,  with  Ac  —  0.8  during  the 
day  and  Ac  —  4.0  ppb  at  night  the  accuracy  with  which 
A  c  is  made  at  the  Pawnee  site  ranges  from  a  =  [(0.07  + 
0.13)  *  100 /Ac]  =  25%  during  the  day  to  a  =  4%  at  night. 


Homogeneity  and  Stationarity 

Field  micrometeorological  measurements  never 
completely  satisfy  the  assumptions  of  horizontal  uni- 
formity and  steady-state  conditions  assumed  in  the 
derivation  of  [7]  and  [8].  Businger  (1986)  has  provided 
a  gross  test  for  these  problems  given  as  follows: 


we; 

w'c' 

dc 

dc 

at 

at 

[A-10] 


where  ht  is  a  scaling  height,  and 
^  =  dz 

at   zj  at 

o 

If  ht  »  z,  the  sampling  height,  then  advection  can  be 
disregarded.  Data  failing  the  test  ht<z  were  eliminated 
from  the  analysis  herein  (see  table  A-2). 
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Sampling  Height 

Based  on  reported  atmospheric  spectra  (Caughey 
1981;  Hojstrup  1981;  Kaimal  et  al.  1972;  Panofsky  and 
Dutton  1984),  measured  nondimensional  frequencies  typi- 
cally range  between:  10"3<  /z/u  <  10.  Given/  =  14.6  Hz 
and  u  values  between  1  m  s"1  and  10  ms"1,  theoretically  a 
sampling  height  of  z  <  10  *  1  /  14.6  —  70  cm  would 
be  an  acceptable  measurement  height  for  the  Pawnee 
site  system.  The  sonic  anemometer  sensing  path  length 
becomes  the  limiting  factor  near  the  surface;  therefore, 
3  meters  was  chosen  as  the  lowest  sampling  height  for 
the  Pawnee  site,  based  on  the  eddy-correlation  experi- 
ence of  Kaimal  (1964)  and  Haugen  et  al.  (1971)  and 
based  on  possible  chemical  reactions  between  03  and 
NOx  below  2  meters  as  discussed  under  Experimental 
Design.  Recently  Volkov  et  al.  (1986)  have  reported 
valid  C02  flux  data  at  1.7  m  height  using  a  20-cm  path 
length  sonic  with  a  4-Hz  sample  rate.  It  is  considered 
that  3  meters  is  well  within  eddy-correlation  technique 
height  limits  reported  in  the  literature. 

Inadequate  Fetch 

Businger  (1986)  cites  the  rule-of-thumb  fetch  re- 
quirement for  micrometeorological  studies  as  equal  to 
lOOz.  The  Pawnee  site  is  relatively  flat,  except  for  the 
north-south  oriented  swale  approximately  400  m  west 
of  the  Pawnee  site  tower.  Analysis  of  measured  zQ  values 
as  functions  of  wind  direction  lead  to  the  initial  conclu- 
sion that  this  swale  is  not  a  problem  (Zeller  et  al.  1989); 
however,  Km  values  reported  herein  are  affected  by 
wind  direction.  Flux  data  gathered  when  the  average 
wind  direction  was  from  the  instrument  shelter  and 
through  the  tower  (210°  ±  20°  azimuth)  were  disre- 
garded for  analysis  because  of  flow  distortion  (Cermak 
and  Horn  1968).  The  3-  and  8-meter  sampling  heights 
seem  to  provide  an  adequate  footprint  to  obtain  repre- 
sentative flux  measurements.  The  surface  canopy  was 
composed  of  sparsely  vegetated  shortgrass  prairie  (mean 
plant  height  -25  cm)  with  leaf  area  indices  between  0.4 
and  0.6  measured  at  the  Pawnee  site  during  a  normal 
growing  season  (Zeller  and  Hazlett  1989). 
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Figure  A-3 — Composite  wind-speed  correction  as  a  function  of 
wind  direction  for  the  Pawnee-site  sonic  anemometers. 

for  the  remaining  wind  directions  involves  the  an- 
emometer itself:  transducer  shadowing.  Sonic  an- 
emometer shadowing  has  only  recently  been  addressed 
in  the  literature  (Businger  1986;  Conklin  et  al.  1988; 
Wyngaard  and  Zhang  1985).  The  shadowing  problem, 
which  results  in  underestimations  in  wind  velocities 
because  of  transducer-generated  flow  vortices,  is  cur- 
rently unresolved  and  the  subject  of  ongoing  research. 

Based  on  the  wind-tunnel  studies  of  Conklin  (1988) 
and  subsequent  tests  that  have  not  yet  been  published 
(Conklin  et  al.  1988;  Baker  1988),  Massman  et  al.  (1990) 
have  constructed  a  composite  flow-distortion  correc- 
tion for  wind-speed  measurements  specific  to  the  ATI, 
Inc.  model  BH-478/B3  sonic  anemometer.  Figure  A-3 
plots  the  composite  correction  multiplier  for  horizontal 
wind  speeds  relative  to  the  Pawnee  site  configuration 
and  true  north  (TN)  wind  direction.  The  vertical  wind 
component  correction  for  shadowing  does  not  seem  to 
be  a  function  of  direction  and  is  taken  as  1.16  based  on 
the  analysis  of  Massman  et  al.  (1990).  A  complete 
verification  of  these  wind  corrections  is  still  wanting. 
However,  based  on  a  comparison  of  the  surface  energy 
balance  components  between  the  Pawnee  site  eddy- 
correlation  measurements  and  measurements  from  the 
nearby  CSU  soil  lysimeter,  the  corrections  can  account 
for  differences  between  the  two  measurement  tech- 
niques (Massman  et  al.  1990). 

Assumptions 


Flow  Distortion  and  Shadowing 

Flow  distortion  refers  to  the  deleterious  effect  of  the 
sensors'  geometric  configuration  and  positioning  on  the 
quantities  being  measured  because  of  local  changes  in 
the  wind-flow  patterns.  Shadowing  refers  to  the  delete- 
rious effect  of  the  flow  around  the  transducers.  As  stated 
previously,  data  associated  with  upwind  directions 
coming  from  the  instrument  shelter  and  support  tower 
towards  the  sensor  arrays  were  discarded.  By  far  the 
most  significant  impact  on  wind-related  measurements 


A  consolidated  list  of  all  the  assumptions  previ- 
ously discussed  or  alluded  to  is  provided  in  table  A-4  for 
completeness.  These  are  assumed  to  be  true  throughout 
the  study  and  in  the  analysis  of  the  data.  The  weakest 
assumptions  are  horizontal  uniformity,  satisfactory 
fetch,  and  nonreactive  trace  gas.  Analyses  of  the  data 
show  some  influence  by  these  factors,  as  discussed  in 
the  main  text.  The  assumption  of  a  logarithmic  profile  in 
the  surface  layer  and  the  resulting  choice  of  4.9  m  as  the 
applicable  height  for  applying  the  results  herein  has  a 
directly  proportional  impact  on  measured 
nondimensional  gradients,  [15a-c]. 
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Table  A-4.  —  Assumptions 


1 .  Stationarity  —  suspect  data  records  are  culled  from  analyses. 

2.  Horizontal  uniformity  —  data  from  wind  directions  030°  -  250° 
were  analyzed  separated  from  250°  -  030°  wind  data. 

3.  Eulerian  point  average  -  ensemble  average  —  a  28.2-minute 
average  is  sufficient  to  measure  averages,  variances,  and 
covariances  with  a  reasonable  degree  of  accuracy. 

4.  No  hot  spots  or  standing  waves. 

5.  Air  density  is  a  constant. 

6.  Coriolis  forces  are  negligible. 

7.  Satisfactory  fetch  —  data  influenced  by  the  swale  west  of  site  is 
analyzed  separately  from  other  data. 

8.  Sensor  separation  is  close  enough  to  provide  valid  covariance 
measurements. 

9.  Ozone  can  be  treated  as  a  nonreactive  trace  gas  between  3  and 
8  meters. 

10.     The  analyses  of  3-  and  8-meter  flux  and  gradient  data  applies  at 
the  geometric  mean  height  of  4.9  meters  —  (dimensionless 
gradients  are  sensitive  to  the  choice  of  z). 


Measurement  Errors 

Based  on  the  uncertainty  interval  results  shown  in 
table  A-5,  eddy  diffusivities  can  be  estimated  with 
slightly  better  precision  than  nondimensional  gradi- 
ents. Sensible  heat  quantities  and  4>h  are  the  best 
measures,  with  10%  and  14%  uncertainty,  respectively, 
followed  by  momentum  with  24%  and  26% ,  and  finally 
ozone  with  40%  and  41%.  The  values  chosen  for  table 
A-5  and  the  uncertainty  calculations  are  representative 
of  daytime  values,  which  result  in  larger  uncertainty 
intervals  compared  to  nighttime  measures.  Daytime 
measurement  uncertainties  for  Kc/Km,  Kh/Kc,  and  Kh/ 

are  ±  41%,  ±  30%,  and  ±  14%,  respectively. 

It  is  impossible  to  absolutely  state  the  error  for  the 
experimentally  determined  values  of  Kq  and  4>q  ex- 
tracted from  the  Pawnee  site  data.  This  is  due  to  the  ever- 
changing  parameters  involved,  and  the  basic  nature  of 

Table  A-5.  —  Measurement  uncertainty 


micrometeorological  field  measurements.  An  "uncer- 
tainty interval"  can,  however,  be  estimated  assuming  a 
linear  error  with  each  independent  variable  (Sandborn 
1981).  For  Kq,  [20] ,  the  uncertainty  interval,  cok,  for  eddy 
diffusivity  can  be  approximated: 
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For  (j)  ,  [15],  the  uncertainty  becomes: 
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Table  A-5  gives  the  numerical  values  used  to  estab- 
lish uncertainty  for  Kq's  and  4>  s  for  £  <  0.  The  values  in 
table  A-5  are  representative  of  the  various  measures  and 
no  more.  They  are  an  attempt  to  objectively  state  the 
uncertainty  in  each  measure.  The  bases  for  CO:^7  in 
table  A-5  are  the  flux  loss  estimates  provided  in  Zeller 
et  al.  (1989);  therefore,  the  relative  w'q'  values  are  1.0. 
For  ozone,  coAq  is  based  on  figure  A-2  and  Aq  on 
average  daytime  Ac  measurements.  For  temperature, 
coAq  is  based  on  the  error  caused  by  solar  heating.  For 
momentum,  coAq-  is  based  on  figure  A-3.  Table  A-5  also 
presents  the  uncertainty  results,  [A-12]  and  [A-13],  in 
terms  of  percent  error.  Uncertainty  for  the  von  Karman 
parameter,  cok,  although  questioned  by  some  authors 
(Businger  1986)  is  not  estimated, 
values  and  percent  uncertainty."1" 


Quantity 

Ozone 

Sensible  Heat 

Momentum 

w  q 

0.3 

(ppb  m  s"1) 

0.05    (°C  m  s"1) 

0.10  (m2s-2) 

0.04 

(PPb) 

0.1  (°C) 

0.2     (m  s"1) 

0.002  (m) 

0.002 

0.002 

coz 

0.003  (m) 

0.003 

0.003 

co'q' 

1.0 

(ppb  m  s1) 

1.0  (°Cms-1) 

1.0  (m2s2) 

Aq 

0.15 

(PPb) 

1.1  (°C) 

0.9     (m  s"1) 

Az 

5 

m) 

5 

5 

z 

8 

(m) 

8 

8 

(VK) 

40% 

10% 

24% 

K/<M 

41% 

14% 

26% 

+  daytime  only 
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APPENDIX  B. 
SPECTRA  AND  COSPECTRA 

Temperature  Spectra 

The  temperature  time-series  sampled  at  40.32  Hz 
using  AIR,  Inc.  FT-1A-T  platinum  resistance  tempera- 
ture sensors  at  heights  8  and  3  meters  between  11:02  and 
11:30  MST  (Mountain  Standard  Time)  on  July  15, 1989, 
provide  an  example  data  set  to  demonstrate  both  the 
response  of  the  sensor  used  and  the  validity  of  sampling 
at  14  Hz.  Fast  Fourier  Transforms  (FFT)  require  data-set 
populations  of  exact  powers  of  2.  In  this  case  the  first 
65,536  (216)  temperature  data  points  were  selected  out 
of  68,235  recorded  for  each  height  for  the  28.2-minute 
sampling  period.  Figure  B-la  is  a  log-log  plot  of  2,048 
frequencies  for  the  8-meter  temperature  data.  The  solid 
curve,  fSe{f)l  d'2 ,  is  the  empirical  spectra  [B-l]  (Kaimal 
et  al,  1972)  multiplied  by  0*  /  0~'2: 


fSB  (/)^[53.4n/&+24n)5/3,         f  <0.15 
02     ~  [24. 4n  /  (1  +12.5n)5/3,  f>0.15 


[B-l] 


where  n  =  fzl  u,  the  nondimensional  frequency. 

The  vertical  alignment  of  Kaimal's  theoretical  spec- 
trum with  that  of  the  data  appears  to  be  an  art.  Several 
authors  "fit"  their  spectra  to  Kaimal's  (Clark  et  al.  1982) 
or  normalize  relative  to  one  spectral  value  in  the  experi- 
ment (Volkov  et  al.  1986).  The  vertical  location  of  data 
depends  on  the  number  of  original  data  points  and 
frequencies  analyzed  as  well  as  the  energy  within  the 
time  series.  The  convention  used  here  is  to  multiply 


each  power  density  by 


q'2(field)/  a2  (spectrum) 


Figure  B-lb  is  the  same  as  figure  B-la  for  the  3-meter 
height.  Both  figures  demonstrate  the  expected  -2/3 
eddy  cascade  in  the  inertial  subrange,  and  both  figures 
show  minimal  aliasing  at  the  high-frequency  end  of  the 
spectra.  To  test  the  routine  sampling  rate,  14  Hz,  every 
third  data  point  was  selected  from  the  8-meter  68,235 
point  temperature  record  introduced  above  to  provide  a 
13.44  Hz  data  set  (i.e.,  40.32/3).  Zeros  were  added  to  the 
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Figure  B-1  a— Temperature  spectra  at  the  Pawnee  site,  8-meter 
data  collected  at  40.32  hz,  1 1 :02-1 1 :30  MST,  July  15, 1989. 
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Figure  B-1b — Temperature  spectra  at  the  Pawnee  site,  3-meter 
data  collected  at  40.32  hz,  11:02-11:30  MST  July  15, 1989. 

resulting  detrended  22,745  data  points  to  provide  a 
32,768  point  record.  Figure  B-lc  is  the  same  as  figure  B- 
la  only  for  13.44  Hz  rather  than  40.32  Hz.  Some  infor- 
mation is  lost  at  the  higher  frequencies  and  the  overall 
spectrum  is  noisier,  even  after  block  averaging;  how- 
ever, the  basic  shape  and  variance  represented  by  the 
lower  frequencies  is  maintained. 

These  typical  spectra  results  and  the  recommenda- 
tion of  at  least  10  Hz  data-acquisition  rate  (Kaimal  1975) 
lead  to  the  conclusion  that  14  Hz  is  a  reasonable  sam- 
pling rate  for  eddy  correlation  and  will  provide  for  the 
capture  of  the  eddys  involved  in  vertical  flux  transfer.  It 
is  important  to  note  that  the  spectra  in  figure  B-la-c  are 
from  a  single  sampling  period  and  are  quite  reasonable 
given  that  [B-l]  is  the  empirical  result  of  many  sample 
period  data  sets. 

Temperature  Cospectra 

The  temperature  cospectra,  i.e.  w'  6'  (n)  heat  flux,  for 
8  and  3  meters  for  23,770  data  points  collected  at 
14.04  Hz  between  13:32  and  14:00  MST  July  7, 1989,  are 
shown  in  figure  B-2 .  The  empirical  cospectrum  of  Kaimal 
et  al.  (1972)  is  depicted  by  the  solid  line: 
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[B-2] 


where  Cw6[f)  =  vertical  velocity  temperature 
cospectrum.  Each  covariance  cospectra  density  was 
adjusted  by  the  multiplier: 


w'0'  (field <re( cospectrum) 


Vertical  Velocity  Spectra 

Velocity  data  were  collected  at  8-  and  3-meter  heights 
using  Applied  Technologies,  Inc.  model  BH-478B/3 
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Figure  B-1c— Temperature  spectra  at  the  Pawnee  site,  8-meter 
data  collected  at  13.44  hz,  11:02-11:30  MST,  July  15, 1989. 
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Figure  B-2a — Temperature  -  vertical-velocity  cospectra  at  the  Paw- 
nee site,  8-meter  data  collected  at  14.04  hertz,  13:32-14:00  MST 
July  7, 1889. 
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Figure  B-2b — Temperature  -  vertical-velocity  cospectra  at  the  Paw- 
nee site,  3-meter  data  collected  at  14.04  hertz,  13:32-14:00  MST 
July  7, 1989. 


sonic  anemometers.  Figure  B-3  shows  the  vertical 
velocity  spectra  for  23,770  data  points  collected  for 
each  height  at  14.04  Hz  between  13:32  and  14:00  MST 
on  July  7,  1989.  Each  set  was  prepared  as  previously 
discussed. 

The  solid  curve  in  figures  B-3a  and  B-3b  is  the 
empirical  curve  for  unstable  and  neutral  conditions 
suggested  by  H0stroup  (1981): 
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[B-3] 


Note  that  the  spectra  in  the  inertial  subrange  of  figure  B- 
3  fall  short  of  the  expected  -2/3  slope.  This  loss  in 
variance  is  due  to  an  analog  filter  installed  by  the  sonic 
manufacturer  to  provide  an  equivalent  10-Hz  signal;  it 
does  not  result  in  any  significant  loss  in  measured  fluxes 
(Kaimal  1989).  The  researcher  has  the  alternative  choice 
of  monitoring  a  digital  signal  directly  from  sonic  sensor. 
The  analog  output  was  used  for  the  flux  measurements 
discussed  herein. 
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Figure  B-3a— Vertical-velocity  spectra  at  the  Pawnee  site,  8-meter 
data  collected  at  14.04  hz,  13:32-14:00  MST,  July  7, 1989. 


0.1 


in 


0.01 


f  =  14.0  hz 

v  =  -3.6 m 
z/L  =  -0.7 
z  =  8  m: 

<5>  .:  Measured  spectra 
 -(Hajstrup  1981) 


0.001 


0.01 


0.1 


10 


100 


fz/u 


Figure  B-3b— -Vertical-velocity  spectra  at  the  Pawnee  site,  3-meter, 
data  collected  at  14.04  hz,  13:32-14:00  MST,  July  7, 1989. 


CAAM  Ozone  Spectra 

Figure  B-4  shows  the  spectra  for  the  CAAM  (ozone) 
data  collected  from  the  8-meter  level  on  July  7,  1989 
between,  13:32  and  14:00  MST.  The  solid  curve  in  figure 
B-4  is  [B-ll  for/Sc(/)/cl.  The  tail  at  the  high  frequency 
end  of  the  spectra  starts  between  1//=  1.5  to  2.0  seconds. 
This  feature  is  aliasing  due  to  the  CAAM  instrument 
response  and  the  effect  of  the  intake  system.  This  aliasing 
feature  was  also  experienced  in  the  3-meter  spectra.  The 
CAAM  signal  aliasing  did  not  seem  to  affect  flux  mea- 
surements, as  presented  next. 
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CAAM  (Ozone)  Cospectra 

The  CAAM  or  ozone  cospectra,  w'c'(n)  ozone  flux, 
for  8  and  3  meters  for  23,770  data  points  collected  at 
14.04  Hz  between  13:32  and  14:00  MST,  July  7, 1989,  is 
presented  in  figure  B-5.  The  solid  curve  is  [B-2]  with 
Cwc(/)  and  c*  replacing  the  values  for  temperatures  as  all 
scalars  are  expected  to  behave  similarly.  The  relative 
location  of  the  empirical  cospectrum  to  the  data 
cospectrum  would  indicate  either  that  the  covariance 
loss  is  greater  with  the  CAAM  instruments  than  it  is 
with  the  platinum  resistance  temperature  sensors  or 
that  there  is  relatively  less  ozone  transport.  Some  loss  is 
expected  and  has  been  noted  before  (Kristensen  and 
Fitzjarrald  1984;  Moore  1986;  Zeller  et  al.  1989). 
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Figure  B-5a — Ozone  -  vertical-velocity  cospectra  at  the  Pawnee  site, 
8-meter  data  collected  at  14.04  hertz,  13:32-14:00  MST  July  7, 
1989. 
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Figure  B-4 — Ozone  spectra  at  the  Pawnee  site,  8-meter  data 
collected  at  14.04  hertz,  13:32-14:00  MST  July  7, 1989. 
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Figure  B-5b — Ozone  -  vertical-velocity  cospectra  at  the  Pawnee  site, 
3-meter  data  collected  at  14.04  hertz,  13:32-14:00  MST  July  7, 
1989. 
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APPENDIX  C 
DATA  PLOTS 
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Figure  C-1 — Unedited  flux  data:  (a)  ozone;  (b)  sensible  heat;  and 

(c)  momentum. 


E 


I 

I 

Hi 

\ 

I 

130                                               150  170 

Julian  date 

190 

- 

i 

V 

ll 

w 

If 

130  150 

Julian  date 

170 

190 

(c) 

I 

m 

liiil'il 

lit 

i 

in 

130  150  170  190 

Julian  date 


Figure  C-2 — Unedited  data:  (a)  ozone  concentration;  (b)  tempera- 
ture; and  (c)  wind  speed. 
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Figure  C-3 — Unedited  data:  (a)  ozone  gradient;  (b)  potential 
temperature  gradient;  and  (c)  horizontal  wind  speed  gradient. 
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Figure  C-4 — Unedited  data:  (a)  total  solar  radiation;  and  (b)  net 

radiation. 
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APPENDIX  D.  TABLE  OF  SYMBOLS 


Defintion 

%  accuracy  expected  for  average  time  T 
pollutant  ozone  concentration 


ozone  concentration  scale  (w'c'/uj 

specific  heat  at  constant  pressure 

atmospheric  cospectrum  of  w  and  q  as  a  function  of  n 

molecular  diffusivity 
nondimensional  frequency  (n  u/z) 
measured  pollutant  (ozone)  flux 

flux  of  q 

acceleration  due  to  gravity 
sensible  heat  flux 

total  cospectral  response  function  for  vertical  wind  w 
and  quantity  q  as  a  function  frequency 

von  Karman  constant 

eddy  diffusivity  of  quantity  q 


-u3T 


Monin-Obukov  length, 


kg  w'O' 

covariance  rotation  matrix 
frequency  (hertz) 

unit  reactor 
sampling  frequency 

specific  density  of  an  arbitrary 
measured  quantity  or  any  quantity 

correlation  coefficient  (p) 

radius 

de  -fdu)2 

gradient  Richardson  number,  g  —  /  6  — 

dz     \  dz  ) 

variance  spectra  density  for  any 
quantity  q  as  a  function  of  n 

time 

temperature 


Unit 

g  m3;  or  ppb 
ppb 

1005  J  kg-1  °K_1 

m  s_1ppb; 

m  s1  °K;  m2  s 2 


cm2  s_1 


gm  m2  s1;  or 
ppb  m  s-1 

gm  m  2  s1;  or 
ppb  m  s 1 


m  s" 


gm  m2  s"1 


cm2  s"1 


m 


5-1 


5-1 


gm  m  3;  ppb; 


or  m  s 


m 


ppb2;  °K2; 
or  m2  s'2 


3C;  or  °K 
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Td(n) 
T» 

T» 
TR(n) 

T(n) 

Tv(n) 

T 

u 

u, 

V 

vd 

w 


w'c' 


WD 

X 

1 

z 
z 

o 

a 

'q 

r 

£ 

v 

e 

e. 

K 
x 

XE 
P 

a2 

c 

a2 

W 


response  function  for  dynamic  frequency  response 
response  function  for  line  averaging 
response  function  for  mismatched  sensors 
response  function  for  high-pass  recursive  filter 
response  function  for  sensor  separation 
response  function  for  volume  averaging 
time  scale  for  averaging  time  to  obtain  a  %  accuracy 
horizontal  wind  speed 


friction  velocity,  V-u'w'  =  ^tq  /  p 
water  vapor  density 
deposition  velocity 
vertical  wind  speed 
ozone  flux 

wind  direction  related  to  true  north 

ith  coordinate  {xl  =  x;x2  =  y,  x3  =  z) 

=  k./kt,  i  =  2,3 

vertical  coordinate 

roughness  length 

K./K 

b  m 

<j)q(measured)/(j>q(empirical) 
adiabatic  lapse  rate 
stability  parameter  (z/L) 
horizontal  wind  rotation  angle 

potential  temperature,  or  verticale  angle  of  mean  wind  vector 


m  s1 
m  s  1 
gm3 
cm  s1 
m  s1 

ppb  m  s"1;  or 
gm  m  2  s"1 

Degree 


m 
cm 


temperature  scale  w'flVU, 
wave  number  in  x  direction 
latent  heat  of  evaporation 
latent  heat  flux 
density  of  air 

temperature  variance 

variance  of  c 

vertical  velocity  variance 


0.0098  °K  m1 

radians 

°K;  or  radians 

°K 

cycles  m1 
2.45xl06J  kg"1 
gm  m2  s1 
gm3 

CK2 

ppb2 


m2  s'2 
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r  momentum  flux;  and  sensor  response  time  gm  m1  s~2;  s 

r  surface  sheer  stress  gm  m 1  s  2 

kz   3  c 

(f>  dimensionless  ozone  gradient,  —  — — 

c*  3z 

kz  d6 

4>h  dimensionless  sensible  heat  gradient  =  —  — 

d>  dimensionless  wind  shear  =  — 

u*  az 

cp  phase  angle  due  to  sensor  mismatch  radians 

vj;q  universal  stability  influence  function 

co  uncertainty  interval  for  q 
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Zeller,  Karl.  1993.  Eddy  diffusivities  for  sensible  heat,  ozone,  and 
momentum  from  eddy  correlation  and  gradient  measurements. 
Res.  Pap.  RM-313.  Fort  Collins,  CO:  U.S.  Department  of  Agricul- 
ture, Forest  Service,  Rocky  Mountain  Forest  and  Range  Experi- 
ment Station.  44  p. 

A  method  developed  to  accurately  measure  trace  gas  mass  (ozone) 
fluxes  and  gradients  from  dry  grassland  canopy  using  surface-layer- 
based  meteorological  towers  to  investigate  the  behavior  of  eddy 
diffusivities  in  the  atmospheric  surface  layer. 

Keywords:  Eddy  correlation,  dry  deposition,  Pawnee  Grasslands,  03, 
micrometeorology,  eddy  correlation,  trace  gas  fluxes,  ozone  deposi- 
tion, eddy  diffusivity,  sensible  heat  flux,  atmospheric  surface  layer, 
Monin-Obukov  Similarity. 
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The  Rocky  Mountain  Station  is  one  of  eight 
regional  experiment  stations,  plus  the  Forest 
Products  Laboratory  and  the  Washington  Office 
Staff,  that  make  up  the  Forest  Service  research 
organization. 

RESEARCH  FOCUS 

Research  programs  at  the  Rocky  Mountain 
Station  are  coordinated  with  area  universities  and 
with  other  institutions.  Many  studies  are 
conducted  on  a  cooperative  basis  to  accelerate 
solutions  to  problems  involving  range,  water, 
wildlife  and  fish  habitat,  human  and  community 
development,  timber,  recreation,  protection,  and 
multiresource  evaluation. 

RESEARCH  LOCATIONS 

Research  Work  Units  of  the  Rocky  Mountain 
Station  are  operated  in  cooperation  with 
universities  in  the  following  cities: 


Albuquerque,  New  Mexico 
Flagstaff,  Arizona 
Fort  Collins,  Colorado* 
Laramie,  Wyoming 
Lincoln,  Nebraska 
Rapid  City,  South  Dakota 
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